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ABSTRACT 

An  experimental  investigation  of  a  wing  boundar\"  layer  subjected  to  periodic  tur- 
bulent flow  at  a  Reynolds  number  of  500,000  was  conducted.  Non-thrusting  turbulence 
pulses  were  generated  at  a  rate  of  50  pulses  per  second  with  a  turbulence  intensity  near 
10%  by  spinning  a  rod  in  the  free  stream  flow  upwind  of  the  test  airfoil. 

Time-var\'ing  velocity  measurements  were  made  at  three  representative  chord  lo- 
cations (laminar,  transitional  turbulent,  fully  turbulent)  at  one  angle  of  attack  using 
single-element  hot-wire  anemometrv'.  The  characteristics  of  the  boundarv'  layer  velocity 
profiles,  turbulence  intensity  profiles,  and  velocity  spectra  and  total  spectral  power  were 
documented. 

The  time-varying  boundary  layer  response  could  be  characterized  by  three  flow  re- 
gimes: undisturbed  flow;  turbulence  pulse;  and  recoverv'  period.  The  boundar\"  layer 
exhibited  a  cyclic  transition  response  varying  from  the  undisturbed  flow  regime  to  the 
turbulence  pulse  regime  back  to  the  undisturbed  flow  regime.  The  turbulence  pulse  was 
found  to  penetrate  into  the  entire  boundary  layer.  The  turbulence  pulse  prompted 
laminarization  of  the  flow  during  the  recover}'  period.  Laminarization  apparently  re- 
sulted from  rapid  acceleration  of  the  near-surface  flow  within  the  boundary  layer  due  to 
turbulence  induced  momentum  transfer  down  through  the  boundar>'  layer  and  local  flow 
acceleration  following  the  velocity  deficit  of  the  turbulence  pulse.  The  eflects  of  the 
periodic  turbulence  pulse  were  most  noticeable  in  the  recover}'  period  at  the 
transitional  turbulent  and  fullv  turbulent  resions  of  the  wine.  Recoverv  was  senerallv 
characterized  by  lowered  turbulence  levels  and  reduction  in  boundarv'  layer  thickness, 
indicating  a  stabilizing  effect  of  the  external  turbulence  pulse. 
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NOMENCLATURE 

c       wing  chord 

F       tunnel  calibration  factor 

f,       cutofT  frequency 

f„       Nyquist  frequency 

fj       sampling  frequency 

K      acceleration  parameter  -tt  ~a~ 
p,       total  pressure 
p       static  pressure 
Ap     differential  static  pressure 
q       test  section  dynamic  pressure 

L-  c 


u'^ 


R,  chord  Reynolds  number 

Ty  turbulence  intensity     ,     - 

U  streamwise  mean  velocity  or  ensemble  average  velocity 

u  streamwise  instantaneous  velocitv 

V  test  section  velocity 

X  streamwise  coordinate 

y  coordinate  direction  normal  to  and  away  from  surface 

S  boundary  layer  thickness 

V  kinematic  viscosity 
p  density 

subscripts 

e  boundar\'  layer  edge  value 

oo  free  stream  value 

superscripts 

'  fluctuating  quantity 

—  time  mean  value 
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I.     INTRODUCTION 

A.     BACKGROUND 

Airfoil  experimental  research  at  low  chord  Reynolds  numbers  has  become  a  major 
field  of  academic,  militar}',  and  civil  interest  in  the  last  few  years  due  to  the  current  de- 
velopment of  low-speed  remotely  piloted  vehicles  (RPVs),  manned  and  unmanned 
ultra-light  aircraft,  human-powered  aircraft,  and  high  altitude  station  keeping  aircraft. 
Low  Reynolds  number  flight  (R,  <  500,000)  is  generally  considered  as  ver\'  low  speed 
flight,  or  night  with  a  short  wing  chord.  Low  Reynolds  number  flight  also  occurs  at  high 
altitudes  where  the  air  density  is  low.  Most  low  Reynolds  number  data  in  the  past  have 
come  from  experiments  on  gliders  and  model  airplanes.  [Ref  1] 

The  aerodynamic  response  of  a  wing  at  low  Reynolds  numbers  is  considerably  dif- 
ferent than  for  Re^'nolds  numbers  greater  than  about  500,000.  At  low  Reynolds  num- 
bers, laniinar  flow  exists  for  a  distance  downstream  from  the  leading  edge.  As  such,  the 
laminar  to  turbulent  transition  occurs  much  further  back  on  the  wing,  often  in  the  ad- 
verse pressure  gradient  region.  The  transition  area  is  often  quite  broad  and  not  well 
defined.  In  addition,  the  transition  process  is  ver\'  sensitive  to  free  stream  disturbances, 
external  turbulence,  and  surface  imperfections. 

Laminar  flow  tends  to  separate  easily  at  low  Reynolds  numbers  especially  in  the 
presence  of  an  adverse  pressure  gradient.  Leading  and  trailing  edge  separation  is  fre- 
quently observed  as  well  as  separation  along  the  airfoil.  Separation  along  the  airfoil  is 
commonly  referred  to  as  the  "lanunar  separation  bubble".  Laminar  separation  bubbles 
greatly  afiect  the  performance  of  the  wing.  The  most  efilcient  low  Reynolds  number 
airfoils  are  designed  to  obtain  maximum  benefit  from  the  separation  bubble;  but  only 
because  the  separation  bubble  cannot  be  completely  eliminated  at  low  Reynolds  num- 
bers. Numerous  factors  also  affect  the  size,  shape,  and  character  of  the  separation 
bubble  including  free  stream  disturbances,  vibration,  surface  condition,  angle  of  attack, 
etc. 

The  critical  nature  of  the  boundary'  layer  and  the  many  factors  that  infiuence  it  have 
prompted  many  recent  investigations  in  the  characteristics  of  airfoils  and  wings  at  low 
Reynolds  numbers.  Although  investigations  under  steady  state  conditions  provide  im- 
portant lift,  drag,  and  airfoil  performance  data,  they  may  not  reliably  predict  the  laminar 
to  turbulent  transition,  leading  and  trailing  edge  separation,  and  formation  and  reat- 


tachment  of  the  separation  bubble  as  these  phenomena  may  be  influenced  to  a  high  de- 
gree by  unsteadiness  in  the  flow.   [Ref.  1] 

Unsteady  low  Reynolds  number  aerodynamic  research  has  recently  been  undertaken 
on  three  fl^onts: 

Airfoil  response  to  high  steady-state  free  stream  turbulence 

Airfoil  response  to  low  to  moderate  frequency  sinusoidal  modulation  of  the  free 
stream 

Airfoil  response  to  transient  free  stream  disturbances 

Meier  and  Kreplin  [Ref  2],  Potter  et  al.  [Ref  3],  and  Kindelspire  [Ref.  4],  among 
others,  have  investigated  airfoil  performance  and  boundary  layer  development  and  be- 
havior in  the  presence  of  high  free  stream  turbulence.  Meier  and  Kreplin  found  that 
both  the  intensity  and  spectral  distribution  (length  scale)  of  the  free  stream  turbulence 
aflected  airfoil  performance.  The  most  significant  changes  were  noted  when  the  length 
scales  were  on  the  order  of  the  boundar}^  layer  thickness.  Potter  et  al.  also  found  that 
the  free  stream  turbulence  intensity  spread  into  the  entire  boundar}"  layer,  but  the  effect 
was  dependent  on  length  scale.  They  also  noted  that  the  free  stream  turbulence  has  to 
penetrate  and  alter  the  velocity  profile  before  producing  any  significant  change  in  the 
boundary  layer  detachment  location.  Kindelspire  found  that  increasing  levels  of  free 
stream  turbulence  caused  correspondingly  earlier  laminar  to  turbulent  transition.  He 
also  conflrmed  that  turbulence  lensth  scales  on  the  order  of  the  boundarv  laver  thickness 
had  the  most  influence  on  boundarv  laver  thickness. 

Brendel  [Ref  5]  performed  a  wind  tunnel  study  with  unsteady  flow  consisting  of  a 
streamwise  sinusoidal  velocity  variation.  Brendel  found  that,  in  general,  the  imposed 
unsteadiness  had  little  eflect  on  the  mean  boundary"  laver  thickness.  The  unsteadiness 
did,  however,  reduce  the  thickness  and  length  of  the  separation  region  associated  with 
the  separation  bubble.  He  also  found  that  lower  frequencies  (large  length  scales)  had 
no  eflect  on  the  overall  flow. 

Howard  [Ref  6]  most  recently  investigated  the  time  varying  characteristics  of  the 
boundary  layer  of  an  airfoil  immersed  in  a  propeller  slipstream  (see  Figure  1).  The 
propeller  wake  subjected  the  airfoil  to  a  transient  free  stream  disturbance.  Howard 
found  that  the  wing  boundar\'  layer  varied  from  laminar  to  turbulent  in  response  to  the 
passage  of  the  propeller  wake.  He  also  found,  in  general,  that  the  effects  of  the  propeller 
slipstream  tended  to  enhance  the  stability  of  the  boundary  layer  and  reduced  the  drag 
coeflicient  in  the  laminar  portion  of  the  slipstream  cycle  below  its  undisturbed  value. 


Laminar 


Turbulent 


Figure  1.      Wing  Immersed  in  Propeller  Slipstream:      [Ref.  6:  p.  17] 


In  particular,  Howard  noted  a  distinct  relaminarization  process  at  work  during  the 
transition  from  turbulent  to  laminar  flow  as  evidenced  by  a  decrease  in  turbulence  and 
an  increase  in  velocity  near  the  surface. 

B.     THESIS  OBJECTIVES 

The  goal  of  this  study  is  to  thoroughly  investigate  the  boundar}'  layer  response  over 
a  wing  to  periodic  free  stream  disturbance  at  the  laminar,  transitional,  and  turbulent 
portions  of  the  wing. 

This  thesis  describes  the  design,  fabrication,  and  testing  of  a  non-thrusting  periodic 
turbulence  generating  device  to  simulate  transient  phenomena  such  as  wind  gust  or  tur- 
bulent wake  disturbances.  The  investigation  consists  of  wind  tunnel  tests  of  an  airfoil 
operating  in  a  periodically  disturbed  free  stream  at  a  chord  Reynolds  number  of  ap- 
proximately 500,000.  iMeasured  data  are  presented  in  the  form  of  non-stationary  mean 
velocity    and    turbulence    time    histories    at    various    heights    above    the    surface;    as 


ensemble-average  velocity  profiles  and  turbulence  profiles;  and  as  velocity  spectra  and 
total  power  plots. 

These  data  will  form  a  base  for  follow-on  research  at  high  angles  of  attack  and  for 
work  within  the  separation  region  of  low  Reynolds  number  airfoils. 


II.     EXPERIMENTAL  EQUIPMENT 

A.     APPARATUS 

The  major  hardware  items  utilized  in  this  investigation  included  the  wind  tunnel,  the 
wing  section  and  wing  mount  assembly,  the  periodic  turbulence  generator,  the  hot-wire 
anemometer,  the  traversing  mechanism,  and  the  data  acquisition  and  processing  system. 
This  chapter  presents  the  specifications,  design,  setup,  and  operating  characteristics  of 
all  equipment  used. 
1.     Wind  Tunnel 

The  Naval  Postgraduate  School  low  speed  horizontal  wind  tunnel  was  employed 
for  this  experiment.    Figure  2  shows  the  tunnel  layout. 

The  tunnel  is  of  single  return  design  and  measures  64  feet  in  length,  25.5  feet  in 
width,  and  supports  a  32-inch  by  45-inch  by  4-foot  test  section.  The  tunnel  was  designed 
by  Aerolab  Development  Company  of  Pasadena,  California  and  has  been  in  continuous 
use  at  the  Xa\'>"  School  since  the  middle  1950's.  [Ref  7] 

The  tunnel  consists  of  five  basic  units:  the  fan  section;  diffusers;  settling  cham- 
ber; contraction  cone;  and  test  section.  The  tunnel  is  constructed  of  wood  throughout 
and  recently  received  a  coat  of  fire  resistant  paint  on  the  interior  surface  to  meet  OSHA 
standards.  The  paint  roughened  the  surface  somewhat;  however,  the  level  of  free-stream 
turbulence  was  not  significantly  worsened  due  to  the  numerous  anti-turbulence  features 
built  into  the  tunnel. 

The  tunnel  is  powered  by  a  100-horsepower  electric  motor  driving  a  three-bladed 
variable-pitch  propeller.  Tunnel  velocity  is  controlled  by  electrically  setting  blade  pitch 
angle  at  the  control  unit  located  adjacent  to  the  test  section.  A  four-speed  Dodge 
transmission  is  installed  between  the  motor  and  the  fan  to  provide  greater  tunnel  velocity 
control  over  the  design  range  of  the  tunnel  (up  to  290  ft/sec).  The  transmission  was  left 
in  3rd  gear  for  all  data  runs  with  a  fan  shaft  speed  of  750  RPM. 

Eight  radial  stator  blades  are  mounted  directly  downstream  of  the  fan  on  the 
drive  shaft  housing.  These  stator  blades  serve  to  straighten  the  fiow  and  remove  the 
swirl  generated  by  the  fan.   These  are  the  first  of  the  installed  anti-turbulence  devices. 

The  diffuser  sections  are  located  on  either  side  of  the  fan  section.  The  gradual 
increase  in  cross-sectional  area  and  the  long  length  of  the  diffuser  sections  gradually  slow 
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the  air  while  introducing  Httle  additional  turbulence  and  keeping  friction  losses  to  a 
minimum. 

The  large  settling  chamber,  located  at  the  left  of  Figure  2,  damps  any  remaining 
flow  unsteadiness.  In  addition,  two  fine  wire  screens  are  mounted  six  inches  apart  in  the 
settling  chamber.  These  serve  to  further  break  up  turbulent  fluctuations  in  the  How 
prior  to  the  air  entering  the  test  section. 

The  contraction  cone  accelerates  the  air  and  delivers  it  to  the  test  section.  The 
large  contraction  ration  of  10:1  helps  assure  low  turbulence  in  the  test  section. 

The  test  section  of  the  low  speed  wind  tunnel  is  equipped  with  a  full  size  re- 
flection plane  mounted  A  inches  above  the  tunnel  floor.  The  reflection  plane  is  perma- 
nently mounted  and  has  an  electrically-driven  turntable  mechanism  in  the  center  (not 
used  in  this  investigation).  The  test  section  is  additionally  equipped  with  full  length 
lighting  located  in  the  corner  fillets  and  has  a  hinged  window  on  either  side.  The  walls 
of  the  test  section  are  slightly  divergent  to  counteract  the  contraction  resulting  from 


boundarv"  laver  srouth.  Breather  slots  are  installed  immediatelv  downstream  of  the  test 
section  to  maintain  the  test  section  at  atmospheric  pressure  and  to  allow  the  tunnel  to 
breath  as  necessary  on  start  up  and  shutdown. 

Two  sets  of  static  pressure  port  rings  are  installed  in  the  tunnel  and  used  to 
measure  the  tunnel  speed.  One  ring  set  is  installed  just  upstream  of  the  contraction  cone 
and  the  other  set  is  installed  at  the  test  section  entrance.  The  pressure  differential  across 
the  two  ring  sets  is  measured  on  a  water  micro-manometer.  This  pressure  differential, 
once  calibrated,  corresponds  to  the  velocity  in  the  test  section.  The  wind  tunnel  velocity 
calibration  procedure  is  covered  in  Chapter  III. 
2.     Wing  Section  and  Wing  Mount  Assembly 

The  wine  section  and  wine  mount  assemblv  are  shown  in  Figure  3.  This  in- 
stallation  was  designed  by  Kindelspire  [Ref  4:  pp.  11-15]  and  a  complete  description  of 
the  apparatus  is  contained  in  his  work.  The  following  paragraphs  highlight  the  pertinent 
aspects  of  the  wing  and  wing  mount  system. 

The  airfoil  tested  was  a  section  of  a  Bell  helicopter  tail  rotor  (model  and  type 
unknown)  having  a  10-inch  chord  and  24-inch  span.  The  airfoil  shape  is  shown  in  Fig- 
ure 4  and  the  section  coordinates  are  given  in  Appendix  A. 

The  airfoil  section  was  stripped  to  bare  metal,  filled,  sanded,  and  finished  with 
high  gloss  black  lacquer  prior  to  installation.  The  lacquer  finish  was  polished  using  fine 
polishing  compound  and  given  a  final  coat  of  automobile  wax  to  bring  it  to  the  highest 
possible  gloss.  A  glossy,  reflective  surface  was  important  as  the  hot-wire  anemometer 
probe  was  positioned  to  within  0.0032  inches  of  the  surface  by  viewing  the  probe  and  its 
reflection  through  a  sighting  level.  The  probe  positioning  procedure  is  discussed  in 
Chapter  III. 

The  wine  mount  assemblv  consists  of  two  8-in  wide  machined  aluminum 
endplates  mounted  vertically  between  the  reflection  plane  and  the  top  of  the  test  section. 
The  wing  section  is  mounted  between  the  endplates  on  rotatable  inserts  that  allow  the 
angle  of  attack  (AOA)  to  be  varied  from  +35  degrees  to  -15  degrees.  Graduated 
markincs  on  the  rotatable  inserts  are  used  for  coarse  settine  of  AOA.  The  acutal  AOA 
is  determined  by  measuring  the  height  of  the  leading  and  trailing  edges  above  the  re- 
flection plane  and  solving  analytically  for  AOA. 

All  data  collected  during  this  experiment  were  taken  1-inch  left  of  mid-span 
(looking  from  the  rear)  due  to  the  presence  of  several  mars  in  the  wing  surface  at  the 
mid-span  line  which  were  caused  by  the  measuring  micrometer  when  the  wing  section 
coordinates  were  documented.    The  1-inch  offset  was  sufllcient  to  remain  free  of  dis- 
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Figure  3.      Wing  and  Wing  Mount  Assembly 


turbances  generated  by  the  surface  imperfectionsl  while  also  remaining  clear  of  any  end 
plate  effects  [Rcf.  4:    pp.  13-14].    The  wing  was  set  at  0  degrees  angle-of-attack  (AOA) 
on  the  inscribed  end  plate  graduations.  This  resulted  in  an  actual  measured  AOA  ofO.TI 
degrees.   All  data  were  collected  at  this  AOA. 
3.     Periodic  Turbulence  Generator 

Periodic  turbulence  pulses  were  generated  by  spinning  a  cylindrical  metal  rod  in 
the  freestream  just  upwind  of  the  wing  section.  The  turbulence  generating  assembly, 
which  resembles  a  propeller,  is  shown  in  Figure  5.  The  mechanism  consists  of  the  base 
support  assembly,  cylindrical  bar  and  hub,  flexible  drive  shaft,  and  motor  and 
timing  position  sensor.  The  support  assembly  and  cylindrical  bar,  hub  were  manufac- 
tured in  the  NPS  Department  of  Aeronautics  and  Astronautics  machine  shop  facilities. 


1  A  general  rule-of-thumb  is  that  disturbance  propagation  is  confmed  to  a   15-degifee  arc 
downstream  of  a  surface  imperfection. 
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Figure  4.      Experimental  Airfoil  Section 


The  support  assembly  was  constructed  out  of  aluminum  stock  and  designed  to 
bolt  to  the  reflection  plane  such  that  the  hub  was  position  2.5  inches  below  and  1  inch 
ahead  of  the  wing  at  center  span.  The  hub  was  placed  well  below  the  wing  in  order  to 
allow  the  disturbances  generated  by  the  hub  section  to  pass  under  the  wing.  The  leading 
edges  of  the  mount  were  also  rounded  and/or  beveled  to  minimize  drag  and  free  stream 
disturbances.  Figure  6  shows  the  wing/wing  mount  and  periodic  turbulence  generator 
installed  in  the  NFS  wind  tunnel. 

The  vertical  strut  of  the  support  assembly  was  bored  to  accept  the  output  end 
of  the  flexible  drive  shaft.  The  drive  shaft  was  held  firmly  in  place  by  a  bracket  mounted 
on  the  backside  of  the  vertical  strut. 

The  bar  and  hub  were  both  machined  out  of  stainless  steel.  The  18-inch  by 
3/8-inch  diameter  bar  was  secured  to  the  hub  by  a  standard  l;8-inch  roll  pin.  The  hub 
was  bored  and  tapped  to  screw  onto  the  threaded  end  of  the  flexible  drive  shaft  extend- 
ing from  the  support  assembly.  The  hub  was  secured  in  place  on  the  threaded  shaft  by 
a  locknut.  In  addition,  the  counterclockwise  direction  of  rotation  was  chosen  such  that 
the  hub  would  tend  to  tighten  during  operation. 

The  38-inch  flexible  drive  shaft  used  to  couple  the  motor  and  the  bar/hub  was 
purchased  at  a  local  hardware  store  as  an  electric  drill  extension.  As  noted  previously, 
the  output  end  had  a  3,'8-inch  threaded  shaft  onto  which  the  hub  assembly  was  mounted. 
The  input  end  of  the  flexible  shaft  was  made  to  fit  a  i;4-inch  drill  chuck.  Use  of  the 
flexible  drive  shaft,  rather  than  direct  drive,  reduced  the  size,  weight,  and  complexity  of 
the  support  assembly  as  it  did  not  have  to  support  the  heavy  drive  motor. 

The  drive  motor  and  position  sensor  are  shown  in  Figure  7.  The  motor  was 
mounted  on  the  tunnel  floor  well  clear  of  the  wing  environment  at  the  downwind  end 
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Figure  5. 


Periodic  Turbulence  Generating  Assembly 


of  the  test  section.  The  motor  is  a  1  4-horsepo\ver  permanent  magnet  variable  speed 
motor.  Motor  speed  could  be  controlled  over  the  range  60-2500  RPM.  An  adapter  with 
a  3  S-inch  threaded  output  shaft  was  manufactured  to  fit  the  1  2-inch  motor  shaft.  A 
I  4-inch  drill  chuck  was  threaded  on  the  adapter  and  used  to  drive  the  flexible  shaft  in 
this  installation.  A  ferro-magnetic  transducer  was  installed  adjacent  to  the  motor  shaft 
to  sense  shaft  speed  and  shaft  position  from  a  toothed  steel  disk  mounted  on  the  motor 
shaft.  The  toothed  steel  disk,  seen  in  Figure  7,  was  modified  by  removing  all  but  two 
of  the  original  30  teeth.  Thus  modified,  an  electrical  pulse  was  generated  twice  each 
revolution  corresponding  to  each  passage  of  the  bar  in  front  of  the  wing  section.  This 
signal  was  used  to  drive  an  electronic  RPM  counter  and  provided  timing  information  for 
the  data  acquisition  system  and  data  analysis  software. 
4,     Hot-\^"ire  Anemometer 

Velocity  and  turbulence  measurements  were  made  using  a  Thermo  System,  Inc 
(TSI)  lFA-100  System  Model  150  Constant  Temperature  Anemometer  and  Model  157 
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Figure  6.      Wing/Wing  Mount  and  Periodic  Turbulence  Generator  Installed  in  NPS 
\^'ind  Tunnel 

Signal  Conditioner.  Figure  S  shows  the  IFA-100  system  installation.  The  lFA-100 
system  is  used  with  touch  panel  entry  of  all  anemometer  and  signal  conditioner  settings 
and  digital  panel  readouts. 

The  Model  150  Anemometer  module  was  operated  in  the  standard  bridge  con- 
figuration with  a  standard  15-ioot  RG-5SA  U  probe  cable.  This  configuration  provided 
extremely  low  noise  output  signals  with  a  typical  frequency  response  ofDC  to  150  KHz 
[Ref  S:  pp.  8-3  to  8-4],  These  parameters  were  important  in  order  to  accurately  measure 
the  low  level,  high  frequency  turbulence  components  of  the  velocity  signal. 

The  hot-wire  signal  was  processed  using  the  Model  157  Signal  Conditioner  to 
improve  the  velocity  resolution  of  the  measurement  system  and  to  provide  low  pass  fil- 
tering. The  unconditioned  hot-wire  output  voltage  over  the  calibration  velocity  range 
varied  from  1.5  volts  at  a  flow  rate  of -4')  ft  sec  to  1.8  volts  at  140  ft  sec.  With  the  given 
analoe  to  digital  voltace  resolution  of  2.44  millivolts,  the  resultin2  vclocitv  resolution 
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Figure  7.      Drive  Motor  and  Position  Sensor 

was.  at  best.  0.8  ft  sec.  By  applying  a  -1.0  volt  offset  and  gain  of  5  to  the  raw  signal  with 
the  signal  conditioner,  the  output  voltage  over  the  calibration  range  was  made  to  vary 
from  2.5  to  4.0  volts  giving  a  velocity  resolution  of  0.16  it  sec--the  iniprovenient  in  re- 
solution being  directly  related  to  the  gain  setting. 

The  low  pass  filter  of  the  Model  157  Signal  Conditioner  was  set  as  close  as 
possible  to  the  Nyquist  cutofl^  frequency.  The  Xyquist  cutolf  frequency,  f^  ,  is  delined 
bv  the  formula: 


f  =  — 


(1) 


where  f^  is  the  digital  sampling  rate  in  Hertz  (Hz).  To  avoid  serious  errors  in  the  data, 
it  is  important  that  the  signal  have  no  significant  components  in  the  frequency  range 
above  f^  prior  to  digital  con\ersion  [Ref  9:  pp.  68-69).  Setting  the  low  pass  filter  to  f^, 
effectively  removes  all  the  unwanted  high  frequency  components. 
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Figure  8.      IFA- 100  Anemometer  System 

A  Dantec  55pll  straight  general  purpose  hot-wire  probe  with  a  5-micron  di- 
ameter platinum  plated  tungsten  wire  was  used  for  all  measurements. 

The  oscilliscope  and  digital  voltmeter,  also  mounted  on  the  IFA-100  equipment 
rack  (Figure  8),  were  used  to  monitor  the  hot-wire  signal  during  each  run.  The  digital 
voltmeter  was  set  to  read  the  mean  DC,  or  steady  state,  component  of  the  signal  and 
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served  to  monitor  the  steady  state  velocity.    The  oscilliscope  displayed  the  AC,  or  tur- 
bulent, component  of  the  signal  and  served  to  monitor  turbulence  in  the  flow. 

The  calibration  procedures  for  the  anemometer  are  covered  in  Chapter  III  of 
this  report. 

5.  Three  Dimensional  Traverser 

Figure  9  shows  the  Velmex  8300  Control/Driver  3-D  traversing  system  em- 
ployed to  position  the  hot-wire  probe  through  the  boundary  layer.  The  traverser  con- 
figuration used  by  this  investigator  was  identical  to  that  utilized  and  described  by 
Kindelspire  (Ref  4:   pp.  21-24]. 

The  traverser  was  mounted  above  the  test  section,  thus  allowing  the  hot-wire 
probe  to  be  positioned  spanwise,  chordwise,  and  vertically  in  relation  to  the  wing  sec- 
tion. Minimum  incremental  movement  in  any  of  the  three  directions  was  0.000125 
inches.  The  traverser  could  be  controlled  either  manually  or  by  software  cominand  from 
a  personal  computer  interface.  An  IBM  PC/AT  was  used  to  control  the  traverser  in  this 
installation.  The  software  routine  used  to  drive  the  traverser  is  explained  in  a  later  sec- 
tion of  this  chapter. 

6.  Analog-to-Digital  Data  System 

Data  collection  was  accomplished  using  a  PC,'AT-type  personal  computer  cou- 
pled to  a  Metrabvte  Model  DAS-16F  analos-to-disital  converter. 

The  PC,' AT- type  personal  computer  used  for  data  collection  and  analysis  was 
configured  with  a  40  megabyte  hard  disk,  two  high  density  disk  drives,  and  operated  at 
a  throughput  speed  of  12  VI Hz.  The  computer  was  additionally  equipped  with  1  mega- 
byte of  random  access  memory.  The  large  storage  capacity  was  necessary  to  handle  the 
volume  of  data  collected  in  this  investigation.  Each  sample  point  yielded  close  to  1 
megabyte  of  raw  data  and  over  150  points  were  sampled. 

The  Metrabyte  DAS-16F  is  a  12-bit  ver\'  high  speed  analog  to  digital  interface 
designed  for  use  with  IBVI  and  IBM  compatible  personal  computers.  The  board  fits 
directly  in  an  expansion  slot  in  the  computer.  The  analog  inputs  were  brought  in 
through  an  external  interface  junction  box.  The  board  operates  at  a  maximum  rate  of 
100,000  samples,  sec  (Hz)  with  up  to  16  single  ended  or  8  difierential  analog  input 
channels.  The  DAS-16F  offers  a  wide  range  of  operating  parameters  including  four 
triggering  options,  three  data  transfer  options,  nine  selectable  input  voltage  ranges,  and 
various  other  input,  output  modes. 

The  DAS-16F  was  configured  to  read  two  difierential  channels  (hot-wire  voltage 
and  timing  pulse)  over  the  input  range  of  0-5  volts.   This  resulted  in  an  analog-to-digital 
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Figure  9.      3-D  Traverser 


resolution  of  2.44  millivolts.    The  hoard  was  operated  at  its  maximum  sample  rate  of 
100,000  Hz  {50,000  llz  per  channel)  during  this  experiment. 

B.     SOFTWARE 

The  computer  played  an  important  role  in  all  phases  of  data  collection  and  analysis 
in  this  investigation.  The  entire  experiment  was  conducted  using  two  PC  AT-type 
computers.  The  software  programs  created  for  each  phase  of  the  investigation  are  de- 
tailed in  this  section.  The  computer  code  for  each  of  the  programs  developed  by  this 
author  is  contained  in  Appendix  C  through  Appendix  F. 
1.     Data  Collection  Software 

The  programs  in  this  section  facilitated  data  collection  and  performed  the  fol- 
lowing control  functions: 

3-D  Traverser  Control  (TR.-WERSE.BAS) 

Ilot-Wirc  Calibration  (HWCAL.BAS) 
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Analog-io-Digital  Board  control  (LABTECH  NOTEBOOK) 
Data  Acquisition  in  the  Boundar\'  Layer  (BLDATA.BAS) 

a.     Program  TRAVERSE.BAS 

The  program  TRAVERSE.BAS  provides  discrete  or  automated  motor 
movement  control  of  the  VELMEX  8300  traverser  through  interactive  keyboard  com- 
mands from  an  IBM  PC/AT  computer.  The  TRAVERSE.BAS  program  written  by 
Kindelspire  [Ref.  4:  pp.  27-30]  was  modified  by  this  author  to  remove  manual  data  entr\- 
routines  and  to  improve  the  accuracy  of  the  motor  stepping  routines. 

Program  TRAVERSE  gives  the  operator  the  option  to  input  individual 
discrete  step  commands  from  the  keyboard  or  to  elect  automated  movement  sequences. 
Discrete  inputs  were  generally  used  to  position  the  hot-wire  in  the  boundary-  layer  prior 
to  each  investigation.  Discrete  movements  could  be  entered  for  any  or  all  of  the  three 
directions  with  the  minimum  movement  being  0.000125  inches  in  any  one  direction. 
Once  the  hot-wire  was  positioned,  the  operator  then  selected  computer  controlled 
movement  to  advance  the  probe  through  the  boundary'  layer  during  the  data  acquisition 
process. 

The  computer  stepped  the  hot-wire  through  the  boundary  layer  according 
to  the  positioning  algorithm  of  Equation  2: 


I       —      I  villi" 


z(j)  =  BLJl-[sm(-| 


■^^(i-1) 


(n-1) 


(2) 


z(j)  =  Position  in  the  boundary  layer 
BL     =  Estimated  boundary'  layer  thickness 
n      =  Ximiber  of  data  points  desired 
j      =  Incremental  data  point 

The  sine  function  based  algorithm  of  Equation  2  is  designed  to  concentrate 
many  smaller  incremental  steps  at  the  bottom  of  the  boundary  layer  and  gradually  in- 
crease the  increment  size  toward  the  top  of  the  layer.  Special  care  was  taken  in  writing 
the  actual  stepping  routine  to  ensure  motor  movement  commands  were  given  in  multi- 
ples of  the  motor  step  interval  (0.000125  inches)  to  preclude  cumulative  roundoff  errors 
throughout  the  boundary  layer. 

b.     Program  HWCAL.BAS 

This  program,  written  by  the  author,  is  used  to  calibrate  the  hot-wire 
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anemometer.    The  program  was  written  in  Microsoft2  QuickBasic  and  compiled  to  gen- 
erate a  stand  alone  executable  program. 

The  hot-wire  is  calibrated  against  a  reference  pitot-static  probe  located  in 
the  test  section  free  stream  adjacent  to  the  hot-wire.  Pitot-static  input  is  in  centimeters 
(cm)  of  water  as  read  on  a  micro-manometer.  These  data  are  manually  entered  on  the 
keyboard  by  the  operator  for  each  calibration  point  while  the  corresponding  hot-wire 
voltase  from  the  IFA-100  Anemometer  is  dieitallv  read  by  the  Metrabvte  DAS-16F 
analog-to-digital  converter  board.  The  operator  also  enters  the  ambient  temperature 
and  pressure  at  the  conclusion  of  each  run.  One  thousand  hot-wire  samples  are  read  and 
averaged  at  each  calibration  point3.  The  hot-wire  was  generally  calibrated  at  1,2,  3,  4, 
6,  8,  and  10  cm  of  water  providing  accurate  calibration  data  over  the  velocity  range 
40-140  ft;sec.  Program  output  consists  of  five  columns  of  data  suitable  for  plotting  us- 
ing most  current  graphing  software  packages: 

Pitot-static  pressure  (cm  of  water) 

Hot-wire  voltage  (volts) 

Free-stream  velocity  (ft/sec) 

(Hot-wire  voltage)  '•'•"  2 

Velocity  **  .5 

The  hot-wire  calibration  procedures  are  covered  in  Chapter  III  of  this  report. 
c.     LA  B  TECH  NO  TEBOOK 

The  commercial  software  package  LABTECH  NOTEBOOK  ''Version  4.3/ 
marketed  by  Laboratory  Technologies  Corporation-*  provided  set  up  and  operating 
control  of  the  DAS-16F  analos-to-disital  board.  L'sina  the  NOTEBOOK  menus,  the 
operator  was  able  to  enter  all  the  operating  parameters  including  sample  rate,  channels 
to  sample,  run  duration,  start  mode,  output  data  file  name  and  output  device,  and  data 
output  format.  NOTEBOOK  allowed  this  investigator  to  take  and  store  up  to  130,000 
data  samples  during  any  one  run  while  operating  at  the  100,000  Hz  rate.  Other  software 
control  programs  appeared  to  be  limited  to  32,000  samples  under  the  same  operating 
conditions.  This  feature  was  paramount,  as  a  minimum  data  capability  of  92,000  dis- 
crete samples  was  required  at  each  point. 


2  Microsoft  is  a  registered  trademark  of  Microsoft  Corporation 

3  The  suggested  analog  to  digital  sample  rate  is  1000  Hz  with  a  run  period  of  1.01  seconds. 

4  LABTECH  is  a  registered  trademark  of  Laboratory  Teclinologies  Corporation. 


17 


d.     Program  BLDATA.BAS 

This  driver  program  controls  the  analog-to-digital  data  acquisition  process. 
Program  BLDATA.BAS  was  written  in  QuickBasic  and  compiled  to  generate  a  stand 
alone  executable  program.  The  user  inputs  the  number  of  data  points  desired  in  the 
boundary  layer  and  a  common  output  file  name.  BLDATA  then  calls  NOTEBOOK 
which  in  turn  initiates  data  collection.  Once  the  data  are  collected  and  the  output  file 
written  by  NOTEBOOK,  BLDATA  renames  the  output  file  giving  it  a  sequentially 
numbered  tag  and  prepares  for  another  data  point  in  the  boundar\'  layer.  At  the  com- 
pletion of  the  run,  the  investigator  has  a  directory'  of  sequentially  named  files  relating  to 
each  data  point. 

2.     Data  Analysis  Software 

The  author  used  the  data  analysis  software  routines  specifically  developed  by 
Johnson  [Ref  10]  to  analyze  the  periodic,  non-steady  velocity  and  turbulence  data  ob- 
tained in  this  investigation.  These  programs  are  not  reproduced  here  due  to  their  length 
and  considerable  documentation.  A  summary,  however,  of  Johnson's  main  programs 
as  well  as  those  wTitten  by  the  author  is  provided. 

a.  Program  CONVERT.FOR 

This  program  was  written  by  the  author  to  convert  the  raw  hot-wire  voltage 
and  position  sensor  voltage  into  velocity  data  and  a  synchronization  pulse.  Input  data 
was  in  ASCII  integer  form  with  the  integer  values  0  to  4095  representing  0  to  5  volts. 
The  bar  sensor  voltage  was  processed  to  return  a  2  at  the  exact  center  of  the  position 
sensor  pulse  and  a  0  elsewhere.  The  hot-wire  voltage  to  velocity  conversion  constants 
are  obtained  from  the  hot-wire  calibration  data  from  program  HWCAL  discussed  pre- 
viously. 

b.  Program  MEANSMU.FOR 

This  program,  written  by  Johnson  [Ref  10:  pp.  158-166],  computes  the  time 
varying  mean  of  a  non-stationary,  periodic  time  series.  The  data  are  first  ensemble  av- 
eraged in  the  time  domain  and  then  transformed  into  the  frequency  domain  via  a  Fast 
Fourier  Transform  (EFT).  The  ensemble  is  then  smoothed  by  applying  a  low  pass  dig- 
ital fiher  to  the  ensemble.  The  ensemble  is  then  transformed  back  into  the  time  domain 
by  an  inverse  EFT  and  written  to  the  output  file.  The  amount  of  smoothing  is  a  func- 
tion of  the  low  pass  filter  cutoff  frequency.  The  lower  the  cutoff  frequency,  the  greater 
the  smoothing  and  vice  versa.  The  amount  of  smoothing  can  be  specified  by  the  user 
through  keyboard  input.  The  methodology  of  selecting  the  low  pass  cutoff  frequency  is 
detailed  in  Chapter  IV, 
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c.  Program  TURBIN.FOR 

Written  by  Johnson  [Ref.  10:  pp.  173-181].  this  program  computes  the  in- 
stantaneous turbulence  intensity  of  a  non-stationary,  periodic  time  series.  It  is  the  tur- 
bulence intensity  counterpart  of  the  non-stationar\'  mean  program  VIEANSMU. 
TL'RBIN  first  calculates  the  time  varying  mean  velocity  of  each  ensemble  by  low  pass 
filtering  in  the  frequency  domain.  TURBIN  then  determines  the  instantaneous  turbulent 
(fluctuating)  component  at  each  point  on  each  ensemble.  The  turbulence  components 
are  then  ensemble  averaged  and  and  written  to  the  output  file.  As  in  MEANSMU.  the 
user  can  specify  the  cutofi"  frequency  of  the  low  pass  filter  used  to  smooth  the  individual 
ensembles. 

d.  Program  PSD.FOR 

This  program,  written  by  Johnson  [Ref  10:  pp.  182-190],  computes  the 
power  spectral  density  of  a  real  discrete  steady  time  series.  This  program  was  for  spec- 
tral analysis  of  the  free  stream  fiow  in  the  test  section. 

e.  Program  ENSPSDA  V.FOR 

This  program,  written  by  Johnson  [Ref.  10:  pp.  191-202],  performs  spectral 
analysis  on  a  short  time  segment  of  an  ensemble  of  non-stationar>",  periodic,  turbulent 
velocity  data.  Whereas  program  PSD  performs  spectral  analysis  on  a  steady  state  signal, 
ENSPSDA  V  investigates  the  time  varying  spectrum  on  selected  portions  of  the  ensemble 
itself.  As  such,  the  author  was  able  to  separately  investigate  the  undisturbed  flow  seg- 
ment, the  turbulence  pulse,  and  the  laminar  flow;recover\'  period  and  characterize  the 
power  spectrum  and  power  content  of  each. 
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III.     EXPERIMENTAL  PROCEDURE 

A.     WIND  TUNNEL  TEST  SECTION  CALIBRATION 

The  purpose  of  wind  tunnel  test  section  calibration  was  to  determine  the  relationship 
between  the  actual  free  stream  dynamic  pressure  (q)  in  the  test  section  and  the  measured 
static  pressure  differential  (Ap)  across  the  contraction  cone.  This  relationship  is  used  to 
indirectly  determine  test  section  free  stream  velocity  by  measuring  Ap. 

For  this  calibration,  a  reference  pitot-static  probe  was  inserted  in  the  center  of  the 
test  section.  A  differential  micro-manometer  was  connected  across  the  total  and  static 
ports  of  the  pitot-static  probe  to  measure  q.  The  relationship  between  the  pitot-static 
probe  differential  pressure,  q,  and  test  section  velocity  is  given  by  Equation  3; 

q  =  (Pi  -  p)  =  -^r  (^) 

p,  =  Total  pressure  measured  by  the  pitot  tube 

p  =  Static  pressure  in  the  test  section 

q  =  Test  section  dynamic  pressure 

p  =  Density  of  air 

v  =  Velocity  in  the  test  section 

The  static  pressure  differential  across  the  contraction  cone,  Ap,  was  measured  by  a 
second  differential  micro-manometer  connected  across  the  two  static  port  ring  assem- 
blies mounted  in  the  wind  tunnel  walls  just  before  and  after  the  contraction  cone.  In 
most  cases  q  is  linearly  related  to  Ap  and  can  be  written: 

q  =  F(Ap)  (4) 

where  F  is  a  constant  called  the  Tunnel  Calibration  Factor  [Ref  7:  pp.  Al-1  to  A 1-5]. 
F  is  determined  experimentally  by  plotting  the  measured  values  of  q  and  Ap.  The  slope 
F  is  given  at  every  point  by: 

q 
F  =   — —   =  Tunnel  Correction  Factor  (5) 

Ap 

Once  F  has  been  determined,  test  section  velocity  can  be  calculated  from  Ap  by: 


/2ApF 
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Calibration  measurements  were  taken  at  numerous  velocities  over  the  range  of  40  - 
140  ft  sec.  The  two  differential  pressures  were  recorded  at  each  calibration  point. 
Figure  10  shows  the  plot  of  q  versus  Ap  for  this  calibration.  The  relationship  is  ob- 
served to  be  linear  and  the  tunnel  correction  factor  is  1.126.  This  value  dilTered  signif- 
icantlv  from  the  calibration  factor  of  1.075  which  has  been  in  common  use  at  NFS  since 
1982.  The  author  speculates  that  the  original  1982  calibration  was  made  prior  to  in- 
stallation of  the  reflection  plane.  Another  possibility  is  that  the  new  calibration  factor 
reflects  the  numerous  small  changes  made  to  the  tunnel  over  the  years. 

B.     HOT-WIRE  ANEMOMETER  CALIBRATION 

The  purpose  of  hot-wire  anemometer  calibration  was  to  determine  the  equation  re- 
lating hot-wire  output  voltage  to  velocity.  This  equation  was  used  in  the  data  conver- 
sion program  CONVDAT.FOR  described  earlier  in  Chapter  III.  The  hot-wire  was 
calibrated  before  and  after  each  data  run.   Hot-wire  calibration  consists  of  three  steps: 

I  FA- 100  system  setup  and  initialization 

Calibration  data  acquisition 

Curve  plotting  and  least  squares  fit 

Initialization  procedures  for  the  I  FA- 100  system  followed  the  steps  for  routine  op- 
eration outlined  in  the  instruction  manual  [Ref  8:  pp.  3-19].  The  setup  parameters  for 
the  I  FA- 100  were  as  follows: 

Bridse  Standard  Bridse  Configuration 

High  Pass  Filter  Not  used 

Low  Pass  Filter  20  KHz 

Voltage  Offset  -1  volt 

Voltage  Gain  5 

As  mentined  in  Chapter  II,  the  low  pass  filter  was  set  to  the  Nyquist  cutoff  frequency 
calculated  by  Equation  1.  The  Nyquist  cutoff  frequency  for  this  investigation  was  25 
KHz.  The  low  pass  filter  could  only  be  set  in  increments  of  10  KHz,  however,  so  the 
filter  was  set  to  20  KHz-the  closest  lower  frequency.  Voltage  offset  and  gain  were  used 
to  improve  the  measurement  resolution  of  the  analog-to-digital  converter  by  expanding 
the  output  voltage  range  of  the  hot-wire  anemometer,  as  discussed  in  Chapter  II. 

The  hot-wire  was  calibrated  against  a  reference  pitot-static  probe  located  near  the 
top  of  the  test  section.  The  hot-wire  was  traversed  to  a  position  adjacent  to  the  pitot- 
static  probe  for  each  calibration. 
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Figure  10.      Wind  Tunnel  Test  Section  Calibration 
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The  hot-wire  was  calibrated  at  or  near  seven  selected  velocities  (1,  2,  3,  4,  6.  8.  and 
10  centimeters  Ap  across  the  contraction  cone)  using  the  software  routine  HWCAL.BAS 
described  in  Chapter  II.  The  calibration  velocities  are  centered  about  the  expected  op- 
erating speed  of  95  ft/sec  and  extended  over  the  range  40  -  140  ft;  sec. 

The  final  calibration  curve  was  obtained  by  plotting  hot-wire  voltage  squared  versus 
the  square  root  of  velocity.  A  nearly  straight  line  relationship  results  for  this  fourth  or- 
der fit.  Most  graphing  routines  will  plot  the  data  and  provide  the  equation  of  the  linear 
least  squares  best  fit  line.   This  linear  equation  is  of  the  form: 

y  =  Ax  +  B  (7) 


where  y  =  ^'velocity   and  x  =  (hot-wire  vohage)^.    Using  the  constants  A  and  B  from 
Equation  7,  the  hot-wire  cahbration  equation  becomes: 

velocity  =  [A(voltage)^  +  B]^  (8) 

This   relationship   can   readily  be   used  in   any   data   conversion   program   such   as 
DATCONV.FOR  (Chapter  II)  to  calculate  velocity. 

C.     HOT-WIRE  POSITIONING 

Great  attention  was  given  to  the  task  of  positioning  the  hot-wire  in  the  boundary 
layer  adjacent  to  the  wing.  The  objective  was  to  accurately  position  the  hot-wire  as  close 
as  possible  to  the  wing  surface  while  not  damaging  either  the  wing  or  the  hot-wire. 

Vaught  and  Fernholz  provide  an  excellent  summary  of  the  various  methods  in  use 
today  [Ref.  1 1].  The  optical  sighting  method  was  used  for  this  investigation.  Kindelspire 
also  used  this  technique  in  prior  work  at  NTS  [Ref  4:  pp.  12-13]. 

Several  modifications  to  Kindelspire's  methods  were  made  by  this  investigator  to 
improve  positioning  accuracy  and  repeatability.  The  techniques  employed  are  described 
below. 

The  height  of  the  hot-wire  above  the  wing  surface  was  measured  by  observing  the 
hot-wire  and  its  reflection  in  the  wing  through  a  surveyor's  sighting  level.  The  sighting 
level  was  positioned  approximately  six  feet  behind  and  to  the  left  of  the  tunnel  test  sec- 
tion. Figure  11  shows  the  sighting  level  in  position  adjacent  to  the  NPS  low  speed  wind 
tunnel. 

As  noted  previously  in  Chapter  II,  the  wing  was  finished  in  high  gloss  black  lacquer 
to  enhance  the  reflected  image  of  the  hot-wire  probe.  The  distance  between  the  hot-wire 
probe  and  its  image  was  determined  by  using  the  horizontal  cross-hair  and  vertical  offset 
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Figure   11.      Surveyor's  Sighting  Level 


vernier  on  the  sighting  level.  The  vertical  vernier  was  adjustable  over  a  range  of  ±  45 
vernier  units  from  the  center  position.  One  vernier  unit  was  determined  to  be  0.0032 
inches  of  vertical  probe  displacement  using  the  traverser,  with  0.000125  inches  stepping 
accuracy,  as  a  calibration  reference.  The  sighting  level  could  be  repeatedly  set  and  read 
to  1  2  vernier  unit  resulting  in  a  positioning  accuracy  of  0.0016  inches. 

The  actual  probe  to  wing  distance  is  1.2  the  distance  observed  between  the  two  im- 
ages. The  hot-wire  images  were  brought  to  within  three  vernier  units  separation  for  data 
collected  at  the  70"  o  chord  location  resulting  in  an  initial  hot-wire  distance  of  0.0048  ± 
O.OOOS  inches.  The  probe  was  brought  to  two  vernier  units  separation  for  all  other  wing 
positions  {50%  and  30- o  chord)  which  equated  to  an  initial  distance  of  0.0032  +  O.OOOS 
inches. 

An  external  flood  light  was  directed  at  the  wing  section  and  the  hot-wire  probe 
during  the  sighting  process  to  enhance  the  image  in  the  wing.  The  poHshed  high  gloss 
black  finish  on  the  wing  provided  a  clear  and  sharp  reficcted  image  comparable  in  qual- 
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Figure  12.      Hot-Wire  Probe  and  Reflection 


ity  to  the  actual  probe  as  shown  in  Figure  12.  No  correction  was  made  for  the  slight 
angle  of  the  wing  surface  due  to  camber  as  this  angle  was  less  than  five  degrees  at  the 
points  examined  resulting  in  a  positioning  error  of  less  than  three  percent. 

D.     ACQUISITION  OF  EXPERIiMENTAL  DATA 

I.     Test  Environipent  Conditions  and  Setup  Parameters 

All  test  were  conducted  at  a  chord  Reynolds  number  of  500,000.    Tunnel  ve- 
locity for  a  siven  Revnolds  number  is  eiven  bv: 


(9) 


R,  =  Chord  Reynolds  Number 

^  =  Viscosity 

p  =  Density 

c  =  Chord  Lensth 
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Solving  Equation  9  using  typical  \PS  ambient  conditions  (temperature  =68  degrees  F, 
pressure  =  29.95  inches  Hg)  where  R,  =  500,000  and  c=  10  inches  yields  a  tunnel  oper- 
ating velocity  of  95  ft  sec.  All  data  runs  were  made  at  this  reference  free  stream  velocity. 
The  setup  parameters  that  remained  constant  throughout  the  investigation  are 
summarized  below: 


Wing  Angle  of  Attack  (AOA) 
Wing  Chord 

Free  Stream  Reference  Velocity 
Chord  Reynolds  Number 
Turbulence  Generator  Motor  Speed 
Frequency  of  Turbulence  Pulses 
Period  Between  Turbulence  Pulses 
Free  Stream  Distance  between  Pulses 
Analog  Channels  Sampled 

Hot- Wire  Voltase 

Timins  Data 
Analos-to-Disital  Rate 
Sample  Period 

Number  of  Pulse  Periods  Sampled 
Number  of  Samples  per  Pulse  period 
Period  Between  Each  Sample 


0.71  degrees 

10  inches 

95  ft;  sec 

500,000 

1,500  RPM 

50  Hz 

20  milliseconds 

22.8  inches 

2  channels 

Channel  1 

Channel  2 

50  KHz  per  Channel 

0.92  seconds 

45 

1000 

20  microseconds 


2.     Equipment  Configuration 

Figure  13  shows  the  equipment  configuration  used  for  this  investigation.  As 
noted  previously,  two  personal  computers  were  used  during  the  data  acquisition  phase. 
One  computer  was  dedicated  to  controlling  the  traverser  mechanism  while  the  other 
managed  data  collection  and  processing.  Two  computers  were  necessarv'  as  the  operat- 
ing programs  and  data  (92,000  samples)  entirely  filled  the  free  memory  of  the  data  ac- 
quisition computer. 

Data  collection  was  quite  routine  once  the  hot-wire  anemometer  was  set  up  and 
calibrated,  and  the  hot-wire  probe  initially  positioned  in  the  boundarv'  layer.  From  that 
point,  the  operator  had  only  to  monitor  the  operating  environment  (tunnel  velocity, 
tunnel  temperature,  turbulence  generator  motor  speed,  and  hot-wire  output),  step  the 
hot-wire  through  the  boundary  layer  using  one  computer,  and  collect  data  with  the  other 
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Figure  13.      Equipment  Configuration 

computer.  Approximately  seven  minutes  were  needed  to  collect  and  store  the  data  at 
any  one  position  in  the  boundary  layer.  While  the  actual  data  sampling  period  was  less 
than  one  second,  the  computer  took  over  5  minutes  to  transfer  the  data  from  random 
access  memory  to  permanent  disk  storage. 

3.     Flow  Characterization 

Data  were  collected  to  characterize  the  flow  in  the  following  regions  of  interest: 

Test  Section  Free  Stream 

Periodic  Turbulence  Pulses 

Boundary  Layer  Response  to  Periodic  Turbulence  Pulses 
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a.  Test  Section  Free  Stream 

Data  were  collected  in  the  free  stream  of  the  test  section  with  all  equipment 
removed  to  determine  the  intensity  and  spectral  content  of  the  flow.  This  step  served 
to  document  the  steady  state  baseline  characteristics  of  the  free  stream  flow  at  95  ft.'sec 
and  also  allowed  this  investigator  to  compare  the  free  stream  turbulence  data  obtained 
in  this  investigation  with  that  obtained  previously  by  other  investigators. 

b.  Periodic  Turbulence  Pulses 

Data  were  collected  with  just  the  periodic  turbulence  generator  installed  and 
operating  in  the  test  section.  The  tunnel  and  the  turbulence  generator  were  set  at  their 
normal  operating  parameters  of  95  ft/sec  and  50  Hz  respectively.  The  hot-wire  was  po- 
sitioned at  the  approximate  vertical  position  of  the  upper  wing  surface  (approximately 
halfway  between  the  hub  and  the  rod  tip).  Measurements  were  taken  four,  six,  and 
eight  inches  downstream  from  the  rotating  rod.  These  three  positions  coincided  with  the 
30%,  50%,  and  70%  chord  locations  on  the  wing.  These  tests  served  to  document  the 
time-varying  velocity,  turbulence,  and  spectral  characteristics  of  the  turbulence  pulse 
generated  by  the  rotating  bar. 

c.  Boundary  Layer  Response  to  Periodic  Turbulence  Pulses 
Measurements  were  made  in  the  boundar\'  layer  at  the  30%,  50%,  and  70% 

chord  locations  on  the  upper  surface  of  the  wing  with  the  periodic  turbulence  generator 
operating  at  50  Hz.  The  wing  AOA  was  measured  at  0.71  degrees  and  the  free  stream 
velocity  was  set  at  95  ft  sec.  Data  were  taken  at  up  to  50  points  over  the  estimated 
boundan.'  layer  thickness  at  each  chord  location.  The  spacing  of  the  sample  points  taken 
at  each  of  the  three  chord  locations  is  given  in  Appendix  B. 

All  data  measurements  were  taken  over  a  four  day  period.  One  additional 
day  of  tunnel  operation  was  needed  following  the  data  collection  runs  in  order  to  inves- 
tigate the  hot-wire  drift  problems  addressed  in  Chapter  IV. 
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IV.     RESULTS 

A.     HOT-WIRE  ANEMOMETER  CALIBRATION 

Hot-wire  calibration  was  conducted  before  and  after  each  boundary  layer  investi- 
gation5  to  confirm  the  stability  of  the  hot-wire  anemometry  system  and  to  ensure  the 
accuracy  of  the  data.  Roane  [Ref  12:  pp.  40-41]  experienced  some  calibration  curve  drift 
in  previous  hot-wire  work  in  the  NFS  low  speed  wind  tunnel.  As  such,  a  minor  amount 
of  drift  was  expected.  This  investigator  planned  to  average  the  initial  and  final  cali- 
bration curves  and  use  the  resulting  average  to  convert  each  day's  data.  The  tabulated 
data  from  the  calibrations  runs  on  each  day  tended  to  confirm  the  expectation  of  only 
slight  drift  and,  therefore,  the  data  acquisition  runs  were  completed  as  planned.  Subse- 
quent detailed  analysis  of  the  calibration  data  for  day  three,  however,  showed  that  the 
initial  and  final  curves  differed  by  approximately  14°  o.  The  drift  for  day  two  was  less 
but  still  significant  at  approximately  8%.  The  drift  for  day  one  could  not  be  determined 
as  the  final  calibration  data  were  lost.  The  results  of  the  calibration  runs  for  day  three, 
presented  in  Figure  14,  show  the  divergence  between  the  initial  and  final  calibration 
data.  The  magnitude  of  the  error  was  too  great  to  resolve  by  averaging.  In  addition, 
the  cause  of  the  drift  and  its  characteristics  were  unknown. 

This  investigator  planned  and  conducted  an  additional  six-hour  tunnel  run  on  the 
followins  dav  in  order  to  document  and  characterize  the  observed  hot-wire  drift.  The 
tunnel  was  run  at  95  ft/sec  with  just  the  pitot-static  and  hot-wire  probes  installed  in  the 
test  section.  Calibrations  were  conducted  at  the  start  and  every  hour  thereafter.  The 
results  of  this  run,  plotted  in  Figure  15,  show  the  shift  of  the  calibration  curve  over  time. 
The  drift  rate  is  greatest  at  the  beginning  and  becomes  very  small  after  six  hours.  The 
effects  of  the  drift  and  the  induced  error  are  greatest  at  the  higher  velocities  due  to  the 
fourth  order  relationship  that  exists  between  hot-wire  voltage  and  velocity. 

The  drift  rate  was  graphically  analyzed  at  a  reference  velocity  of  100  ft/ sec  and  ex- 
panded by  extrapolation  to  cover  an  eight-hour  period.  Figure  16  shows  the  change  in 
calibration  of  the  hot-wire  over  eight  hours.  Notice  that  nearly  50%  of  the  drift  occurs 
in  the  first  two  hours  of  operation. 


5  The  final  calibration  data  for  day  one  was  not  saved  due  to  operator  error. 
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Figure  14.      Hot-Wire  Calibration  Results,  Day  Three 
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Figure  15.      Hot-Wire  Drift  Analysis 
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Figure  16.      Hot-Wire  Calibration  Drift  Over  Eight  Hour  Period 
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The  author  applied  the  drift  vs.  time  relationship  of  Figure  16  to  each  day's  cali- 
bration data  in  order  to  generate  a  set  of  time  var>"ing  calibration  equations.  The 
equations  were  staggered  such  that  the  maximum  velocity  diflerence  between  sequential 
calibration  curves  was  one  ft.'sec  or  less.  These  corrections  assumed  that  the  drift  rate 
determined  during  the  test  run  was  representative  of  the  drift  experienced  on  each  day 
and  that  the  run  duration  for  each  days  was  eight  hours.  Corrections  were  applied  only 
to  the  data  taken  on  day  two  (50%  chord)  and  day  three  (30%  chord). 

This  investigator  believes  the  hot-wire  drift  is  a  result  of  environmental  contam- 
ination  of  the  hot-wire  probe.  Hot-wire  contamination  appears  to  be  a  well  known 
phenomenon  among  experienced  researchers,  but  a  subject  about  which  little  is  printed 
in  current  hot-wire  operating  manuals  or  in  fluid  flow  measurement  texts.  An  older 
publication  by  TSI  has  this  to  say  about  hot-wire  stability: 

For  a  mean  flow  measurement  to  be  useful  it  must  maintain  its  calibration  over 
long  time  periods.  Flot  wires  are  notoriously  troublesome  in  this  regard  due  to  their 
small  size  making  them  susceptible  to  picking  up  small  contaminants  ....  The 
larger  hot  film  sensor  has  been  a  great  improvement  in  this  regard.  The  elTects  of 
surface  contamination  tend  to  be  a  function  of  sensor  size  and  a  ceramic  coating  on 
the  sensor  surface  prevents  corrosion  and  erosion.  Hence,  the  hot  film  is  far  supe- 
rior for  mean  measurements.  [Ref  13:  p.  91] 

Figure  17,  reprinted  from  the  above  TSI  publication,  shows  the  typical  change  in 
calibration  of  a  hot-wire  sensor  and  a  hot-film  sensor  vs.  time  w'hen  subjected  to  an  ar- 
tificially dirty  flow  stream.  The  drift  vs.  time  curve  for  this  investigation  (Figure  16)  is 
remarkably  similar  to  the  TSI  drift  curve.  With  this  knowledge,  it  was  not  difficult  to 
pinpoint  the  source  of  contamination  to  the  wind  tunnel  itself 

For  the  month  prior  to  this  investigation,  the  tunnel  was  used  for  a  series  of  flow 
visualization  studies  using  an  oil-based  smoke  as  a  visualization  medium.  This  type  of 
medium  tends  to  produce  a  certain  amount  of  contaminants  which,  due  to  the  closed 
circuit  design  of  the  tunnel,  were  trapped  in  the  tunnel  environment.  This  oil  residue  is 
the  likely  contaminant  and  cause  of  the  observed  drift. 

Little  can  be  done  about  the  tunnel  contamination,  but  one  likely  solution  to  the 
drift  problem  is  the  use  of  a  hot-film  sensor  instead  of  a  hot-wire  probe.  As  noted  pre- 
viously, the  hot-film  is  virtually  impervious  to  particulate  contamination.  A  second 
solution  would  be  to  pre-run  the  tunnel  prior  to  taking  data  for  six  to  eight  hours  to 
allow  the  hot-wire  to  stabilize. 
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Figure  17.      Change  in  Calibration  of  Hot-Wire  vs.  Hot-Film  When  Subjected  to  a 
Dirty  Flow  Stream:     [Ref.  13:  p.  92] 

B.     ANALYSIS  OF  TIME  VARYING  DATA 

Ensemble  averaging  methods  were  employed  to  characterize  the  non-stationary,  pe- 
riodic data  acquired  during  this  investigation.  The  periodic  nature  of  the  data  made  it 
possible  to  collect  a  series  of  repetitive  time  histories,  or  pulse  passage  cycles,  each  as- 
sociated with  the  passage  of  one  turbulence  pulse.  This  series  of  time  histories  is  called 
an  ensemble.  A  typical  ensemble  consisted  of  a  series  of  45  pulse  passage  cycles.  Ihe 
number  of  passage  cycles  in  one  ensemble  was  essentially  hmited  by  available  computer 
memory.     For  this  investigation,  45  cycles  was  the  upper  limit.    Howard  [Ref  6:  pp. 
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114-115]  showed  that  20  or  more  pulses  were  needed  in  each  ensemble  to  accurately 
portray  the  actual  mean  values,  substantiating  that  45  pulses  were  sulTicient  in  this  case. 

1.  Ensemble  Characteristics 

Each  of  the  45  pulse  passage  cycles  making  up  the  typical  ensemble  for  any  one 
data  point  is  slightly  different  due  to  the  random  nature  of  the  turbulent  components 
of  the  signal.  Figure  18  and  Figure  19  show  two  sequential  pulse  passage  cycles  (in- 
stantaneous velocity  vs.  sample  number)  from  a  typical  ensemble  (Test  Case  A)  taken 
near  the  wing  surface  at  10%  chord.  These  data  were  collected  using  a  single  hot-wire 
sensor.  The  velocity  data  were  sampled,  or  digitized,  at  50  KHz  with  the  turbulence 
generator  operating  at  50  Hz.  This  digitizing  rate  results  in  approximately  1000  data 
points  per  passage  cycle.  The  period  of  each  passage  cycle  being  20  milliseconds,  each 
sample  represents  20  microseconds  in  time.  Time  is  not  recorded  expHcitly  but  can  be 
calculated  from  the  digitizing  rate  and  sample  number.  The  data  are  plotted  with 
straight  line  segments  between  each  of  the  sample  points. 

Each  passage  cycle  is  typically  made  up  of  three  representative  flow  regimes  as 
depicted  in  Figure  18: 
Undisturbed  Flow- 
Turbulence  Pulse 
Recovery  Period 

The  specific  characteristics  of  these  regimes  are  ver\'  dependent  on  chord  location  and 
position  in  the  boundary'  layer;  but  the  following  general  characteristics  hold  true  for 
each.  The  undisturbed  flow  segment  is  representative  of  the  natural  steady-state  flow  for 
the  wing  when  operated  in  an  undisturbed  free  stream.  Howard  [Ref.  6:  p.  81]  confirmed 
this  relationship  experimentally  by  comparing  the  drag  of  the  undisturbed  segment  with 
that  for  natural  flow  and  finding  the  two  equal.  The  turbulence  pulse  is  clearly  identified 
as  a  highly  turbulent  velocity  deficit.  The  recovery  period  is  broadly  defined  as  the  pe- 
riod from  the  end  of  the  turbulence  pulse  to  the  point  where  the  flow  returns  to  the  un- 
disturbed state. 

2.  Mean  Velocity  Determination 

The  author  employed  program  MEANSMU.FOR  described  in  Chapter  II  to 
extract  the  slowly  varying  mean  velocity  from  the  ensemble  of  instantaneous  velocities. 
In  principle,  .VIEANSMU  filters  out  the  high  frequency  turbulence  components  while 
leaving  the  low  frequency  variations  of  the  signal  mean  intact.   The  amount  of  smooth- 
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Figure  18.      Pulse  Cycle  Number  One  for  Test  Case  A:    (Fs=  50  KIlz) 
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Figure  19.      Pulse  Cycle  Number  Two  for  Test  Case  A:    (Fs=  50  KI  Iz) 
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ing  is  determined  by  the  user  selectable  cutoff  frequency  of  the  digital  low  pass  filter  in 
program  MEANSMU. 

Program  MEANS  ML'  extracts  the  mean  through  a  two  step  process.  First,  the 
45  time  histories  are  ensemble  averaged  to  generate  a  single  "representative"  time  his- 
tory. Ensemble  averaging  consists  of  simply  stacking  the  45  time  histories  one  on  an- 
other and  taking  the  average  of  the  45  values  at  each  sample  point  in  the  time  history. 
Figure  20  presents  the  ensemble  average  of  the  45  pulse  passage  cycles  for  Test  Case 
A.  Note  that  the  large  turbulent  fluctuations  observed  in  Figure  18  and  Figure  19  have 
been  reduced  but  not  entirely  eliminated.  The  three  flow  regimes  are  also  quite  pro- 
nounced and  easily  identifiable. 

Figure  21  through  Figure  24  show  the  application  of  digital  low  pass  filtering 
to  the  ensemble  average  of  Figure  20  with  low  pass  filter  cutoff  frequencies  (fj  of  1082 
Hz,  541  Hz,  360  Hz,  and  267  Hz  respectively.  Note  that  the  amount  of  smoothing  is 
directly  related  to  f,. ;  the  lower  f,.  the  greater  the  smoothing.  360  Hz  was  choosen  as  the 
optimal  cutoff  frequency  based  on  the  findings  of  Landrum  and  Macha  [Ref  14]  in  their 
analysis  of  the  spectra  documenting  the  growth  of  turbulent  disturbances  through  the 
transition  process.  They  found  that  the  free  stream  energy  is  contained  below  200  Hz 
while  the  turbulence  enerev  is  generallv  contained  above  550  Hz  in  the  earlv  transition 
stages.  The  cutoff  frequency  was  choosen  to  lie  between  these  two  frequencies. 
Coincidentially,  the  f^=  360Hz  curve  is  relatively  smooth  and  flat  throughout  the  undis- 
turbed flow  regime  and  displays  good  peak  response  over  both  the  turbulence  pulse  and 
recovery^  period  regimes. 

3.     Turbulence  Level  Determination 

Turbulence  values  of  the  time-varying,  periodic  signals  were  calculated  using  the 
program  TLRBIN.FOR  described  in  Chapter  II.  TLRBIN  accomplishes  this  by  calcu- 
lating the  instantaneous  turbulence  at  each  point  in  each  of  the  45  time  histories  and 
then  ensemble  averaging  the  instantaneous  turbulence  magnitudes.  The  result  is  a 
time-varying  representation  of  the  mean  turbulence  at  each  sample  point  in  the  pulse 
period.  It  is  important  to  note  that  TLRBIN  also  employs  user  selectable  digital  low- 
pass  filtering  in  its  routines.  The  author  used  360  Hz  in  this  routine  to  maintain  conti- 
nuity among  the  various  output  data. 

Figure  25  shows  the  turbulence  trace  for  Test  Case  A.  Note  that  the  represen- 
tative flow  regimes  (undisturbed  flow,  turbulence  pulse,  and  recovery  period)  are  as 
clearly  identifiable  here  as  in  the  mean  velocity  traces  in  the  previous  section. 
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Figure  20.      Ensemble  Average  for  Test  Case  A:     (Fs  =  50  KHz) 
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Figure  21.      Smoothed  Ensemble  Average  for  Test  Case  A,  Fc=  1082  Hz:     (Fs  =  50 
KHz) 
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Figure  22.      Smoothed  Ensemble  Average  for  Test  Case  A,  Fc  =  541  Hz:    (1-8  =  50 
KHz) 
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Figure  23.      Smoothed  Ensemble  Average  for  Test  Case  A,  Fc  =  360  Hz:    (Fs  =  50 

KHz) 
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Figure  24.      Smoothed  Ensemble  Average  for  Test  Case  A,  Fc  =  267  Hz:    (15  =  50 

KHz) 
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Figure  25.      Smoothed  Mean  Turbulence  for  Test  Case  A:    (Fc=360  Hz,  Fs  =  50 
KHz) 
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4.     Spectral  Analysis 

Spectral  analysis  of  the  time-van'ing,  periodic  signals  was  conducted  using  pro- 
gram ENSPSDAV.FOR  described  in  Chapter  II.  EXSPSDAV  performs  spectral  analy- 
sis on  a  user  selectable  256  sample-wide  window  of  the  pulse  passage  cycle  time  histor>'. 
In  theor\",  spectral  analysis  can  only  be  accomplished  on  a  discrete  (digitized)  time  series 
of  infinite  length.  In  practice,  the  window  width,  or  number  of  samples,  can  be  reduced 
through  the  application  of  special  techniques  to  the  Fast  Fourier  Transform--the  com- 
promise being  loss  of  frequency  resolution  with  reduction  in  window  width.  The  256 
sample-wide  window  used  in  this  investigation  was  selected  as  the  best  balance  of  these 
factors.  Each  flow  regime  fits  into  the  window  and  the  resulting  frequency  resolution  is 
190  Hz.  The  reader  must  keep  in  mind  that  the  output  data  from  EXSPSDAV  repres- 
ents the  average  taken  over  the  entire  window.  Thus,  dynamic  changes  that  occur  in  the 
flow  regime  over  the  window  period  will  be  lost  in  the  resulting  average. 

Figure  26  and  Figure  27  show  the  output  from  program  EXSPSDAV  for  the 
turbulence  pulse  regime  of  Test  Case  A.  Figure  26  displays  the  velocity  spectrum  of  the 
signal.  The  vertical  scale  is  a  dB  measurement  of  the  relative  turbulence  energy  with  0 
dB  representing  the  highest  magnitude  component  of  the  signal.  The  horizontal  scale  is 
frequency  in  Hz  plotted  out  to  the  Xyquist  frequency.  Each  10  dB  reduction  in  the  ve- 
locity spectrum  represents  a  10  fold  reduction  in  power  level.  For  example,  the  power 
level  at  4.500  Hz  (-20  dB)  is  1,100  the  level  of  the  highest  magnitude  component.  The 
mean  has  been  removed  from  the  signal6  so  that  only  the  spectrum  above  0  Hz  is  ana- 
lyzed. It  is  important  to  note  that  the  spectrum  display  provides  no  information  about 
the  absolute  power  contained  in  the  spectrum,  but  only  displays  power  levels  relative  to 
the  highest  magnitude  in  the  spectrum.  Figure  27  displays  the  total  (cumulative)  power 
distribution  plotted  with  regard  to  frequency.  The  total  power  curve  asymptotically 
approaches  the  value  of  the  total  power  in  the  spectrum.  The  shape  of  this  curve  is  in- 
dicative of  the  power  distribution  within  the  spectrum.  The  sharp  knee  observed  in  this 
plot  indicates  that  the  low  frequency  components  of  the  signal  contain  most  of  the 
spectral  power  with  ver\'  little  power  in  the  higher  frequency  components. 

C.     FREE  STREAM  ANALYSIS 

The  results  of  the  free  stream  analysis  are  presented  for  flow  at  95  ft;  sec  with  all 
equipment  removed  from  the  test  section.    Figure  28  is  a  plot  of  velocity  response  over 


6  The  mean  is  removed  for  all  the  velocity  spectrum  and  total  power  distribution  plots  in  this 
report. 
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Figure  26.      Velocity  Spectrum  for  Test  Case  A,  Turbulence  Pulse 

1000  samples.  The  signal  is  flat  with  a  high  frequency  turbulent  component  of  small 
amplitude  superimposed  on  the  mean.  The  free  stream  turbulence  was  calculated  at 
0.182%  by  the  data  analysis  methods  of  this  investigation.  Ihis  figure  compares  favor- 
ably with  the  free  stream  turbulence  intensity  of  0.23%  observed  by  Roane  [Ref  12:  p. 
41]  at  a  velocity  of  116.4  ft/sec.  The  slightly  lower  value  of  free  stream  turbulence  ob- 
served here  is  not  unexpected  due  to  the  slower  free  stream  velocity. 

Figure  29  shows  the  velocity  spectrum  of  the  free  stream.  The  curve  shows  a  sharp 
drop  in  intensity  with  increasing  frequency  indicating  that  the  major  frequency  compo- 
nents are  at  the  lower  frequencies. 
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Figure  27.      Total  Power  Distribution  lor  Test  Case  A,  Turbulence  Pulse 

The  total  power  distribution  of  the  free  stream  is  depicted  in  Figure  30.  The  rapid 
rise  and  sharp  knee  of  this  curve  confirm  the  concentration  of  power  in  the  lower  fre- 
quencies. The  frequency  components  above  2  KHz  provide  Httle  additional  power  as 
shown  by  the  flatness  of  the  curve  above  this  point. 

D.    TURBULENCE  PULSE  ANALYSIS 

The  turbulence  pulse  was  examined  at  three  positions  (4",  6",  and  8 ")  downstream 
of  the  spinning  rod  with  the  wing  and  wingmount  assembly  removed.  The  smoothed 
mean  velocity  time  histories  for  the  three  positions  are  shown  in  Figure  31.  The  plots 
show  a  nearly  symmetric  velocity  deficit  trace  at  all  three  positions.  Pulse  duration  is 
roughly  250  samples  or  5  milliseconds  and  the  deficit  maxima  are  constant  at  -20  ft/sec. 
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Figure  28.      Test  Section  Free  Stream  Flo>v:    (Fs=  50  KHz) 
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Figure  29.      Free  Stream  Velocity  Spectrum 

The  width  of  the  delicit  appears  to  increase  slightly  with  increasing  distance  downstream. 
The  velocity  deficit  trace  observed  for  the  spinning  rod  is  very  similar  to  that  generated 
by  a  stationary  cylinder  in  steady  cross  flow.  This  similarity  shows  that  spinning  the  rod 
docs  not  significantly  alter  the  classical  response  of  the  cylinder. 

Figure  32  is  the  plot  of  the  smoothed  mean  turbulence  traces  of  the  turbulence  pulse 
at  the  same  positions.  The  pulse  is  intensely  turbulent  with  a  peak  value  of  10,5°o.  The 
turbulence  traces  coincide  with  the  velocity  traces;  however,  the  turbulence  traces  be- 
come somewhat  non-symmetric  downstream.  In  addition,  the  turbulence  appears  to 
anticipate  or  lead  the  velocity  deficit  by  up  to  60  samples  or  1.2  milliseconds.  .Most 
notably,  the  turbulence  magnitude  decreases  quite  rapidly  with  distance  downstream. 
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Figure  30.      Free  Stream  Total  Po>ver  Distribution 


This  decay  is  not  unexpected  and  is  generally  attributed  to  dissipation  of  the  high  fre- 
quency turbulence  components  due  to  viscosity  [Ref.  15:  p.  48]. 

Figure  33  shows  the  velocity  spectra  of  the  turbulence  pulse  regime  and  the  undis- 
turbed flow  regime  taken  at  4"  downstream  from  the  spinning  rod.  The  turbulence  pulse 
trace  decays  more  rapidly  and  to  a  lower  level  than  the  undisturbed  flow  trace.  Addi- 
tionally, the  undisturbed  flow  trace  levels  out  at  approximately  -20  dB  for  frequencies 
above  5  KHz  while  the  turbulence  pulse  trace  bottoms  out  at  -42  dB.  These  plots  show 
the  turbulence  pulse  to  be  a  low  frequency  dominated  event  while  the  undisturbed  flow 
regime  has  components  at  much  higher  frequencies. 
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Figure  31.      Turbulence  Pulse  Mean  Velocity:    (Fs  =  50KHz) 


51 


6  6  6 

td  (d  cd 

0)  Q)  (U 

u,  ^  u 

^-»  -►J  -*-» 
www 

0  ^  ^ 

^  ^  ^ 

o  o  o 

Q  Q  Q 


(-juaojad)    nj^ 


Figure  32. 


Turbulence  Pulse  Mean  Turbulence:    (Fs  =  50  KHz) 
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Figure  33.      Turbulence  Pulse  Velocity  Spectrum 

The  total  power  curves  for  the  two  regimes  are  plotted  in  Figure  34.  1  he  intensity 
of  the  turbulence  pulse  is  clearly  shown  here;  the  pulse  power  is  4  orders  of  magnitude 
greater  than  the  undisturbed  How.  The  rapid  rise  and  sharp  knee  of  the  pulse  trace  also 
confirms  the  intense,  low  frequency  nature  of  the  turbulence  pulse.  Very  little  power  is 
added  by  the  frequency  components  above  1  KHz. 

E.     BOUNDARY  LAYER  ANALYSIS 

The  response  of  the  boundar>'  layer  to  periodic  turbulence  was  investigated  at  three 
locations  on  the  upper  surface  of  the  wing:  30%  chord,  50%  chord,  and  lOVo  chord. 
These  three  locations  were  chosen  as  the  natural  steady-state  (undisturbed)  flow  was 
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Figure  34.      Turbulence  Pulse  Total  PoMer  Distribution 
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laminar  at  30%  chord,  transitional/turbulent  at  50°/o  chord,  and  fully  turbulent  at  70% 
chord. 

1.     30%  Chord 

The  ensemble-average  velocity  and  turbulence  profiles  are  shown  in  Figure  35. 
These  plots  consist  of  a  series  of  time-sequenced  profiles,  or  "snapshots",  taken  at  several 
points  in  the  time  history.  These  profiles  depict  the  dynamic  characteristics  of  the 
boundary  layer  throughout  one  pulse  passage  cycle. 

The  profiles  show  that  the  flow  at  30%  chord  transitions  from  undisturbed 
laminar  flow  with  very  little  turbulence  {sample  100);  to  highly  turbulent  flow  with  a 
nearly  constant  level  of  turbulence  at  the  turbulence  pulse  center  (sample  275);  and  back 
again  to  laminar  flow  (sample  550).  It  is  interesting  to  note  how  quickly  and  uniformly 
the  boundary  layer  responds  to  the  turbulence  pulse.  The  turbulence  profiles  show  a 
nearly  uniform  turbulence  increase  and  corresponding  decrease  with  the  passing  of  the 
turbulence  pulse.  At  the  same  time,  the  velocity  profiles  transition  quickly  from  laminar 
to  turbulent  and  back  to  laminar  in  step  with  the  turbulence  pulse.  The  velocity  increase 
observed  near  the  surface  of  profile  350  is  unusual.  The  velocity  increase  occurs  just 
after  the  turbulence  pulse  and  is  still  evident  in  the  following  profile  (sample  400).  The 
possible  reasons  for  this  behavior  will  be  examined  later  in  this  section. 

The  smoothed  mean  velocity  time  histories  and  the  smoothed  mean  turbulence 
time  histories  of  the  pulse  passage  cycle  are  shown  in  Figure  36  and  Figure  37  respec- 
tively. These  traces  represent  the  time  histories  at  increasing  vertical  distances  (y)  above 
the  wing  surface.  The  bottom  trace  shows  the  flow  history  taken  closest  to  the  wing 
surface  (near-surface  region)  while  the  top  trace  depicts  external  How  outside  of  the 
boundar}'  layer. 

The  velocity  and  turbulence  traces  in  the  external  flow  region  are  typical  of  the 
response  observed  downstream  of  the  periodic  turbulence  generator  (Figure  31  and 
Figure  32);  i.e.,  a  symmetric,  bell-shaped  pulse  followed  by  an  immediate  return  to  the 
free  stream  or  undisturbed  state. 

The  velocity  time  history  traces  show  a  velocity  deficit  of  11  ft/sec  in  the  ex- 
ternal flow  due  to  the  turbulence  pulse.  This  deficit  increases  as  the  surface  is  ap- 
proached until  very  close  to  the  wing.  Quite  near  the  surface,  the  velocity  deficit 
decreases  significantly  and  is  followed  by  a  sharp  velocity  increase  just  after  the  turbu- 
lence pulse.  This  velocity  increase  remains  for  some  time  and  is  slow  to  dissipate.  The 
turbulence  traces,  on  the  other  hand,  show  that  the  turbulence  intensity  of  the  pulse  is 
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Figure  35.      Ensemble  Average  Profiles,  30%  Chord 
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Figure  36.        Mean    Velocity    Time    Histories    in    the    Boundary    Layer,    30% 
Chord:     (Fs=50KHz) 


57 


10-1 


13 


External    Flow 
y   =    0.1731    in. 


y   =    0.0944    in. 


y   =    0.0372   in. 


y   =    0.0184    in. 


y   =    0.0071    in. 


Near   Surface 
y   =    0.0032   in. 


-1 1 1 1 — I — I — I 1 1 — I 1 1 \ 1 1 1 1 1 1 1 

0     100    200    300    400    500    600    700    800    900    1000 


Sample    Number 


Figure  37.       Mean    Turbulence    Time    Histories   in   the   Boundary   Layer,   30% 
chord:     (Fs=50KHz) 
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felt  uniformly  throughout  the  boundar}'  layer  with  only  a  slight  reduction  in  magnitude 
and  stretching  of  the  pulse  very  near  the  surface. 

Figure  38  through  Figure  42  show  the  velocity  spectrum  and  total  power  dis- 
tribution plots  for  the  three  major  flow  regimes  (undisturbed  flow,  turbulence  pulse,  and 
recovery  period)  at  points  in  the  boundary'  layer.  Recall  that  the  velocity  spectrum  plot 
display  is  relative  to  the  highest  component  in  the  spectrum  while  the  total  power  dis- 
tribution plot  depicts  the  cumulative  actual  power  with  regard  to  frequency.  Recall  also 
that  the  mean  has  been  removed  from  the  spectral  analysis.  Hence,  power  changes  in 
the  mean  due  to  mean  velocity  increases  or  deficits  are  not  reflected  in  the  spectral 
curves.  One  must  view  the  mean  velocity  and  turbulence  time  histories  along  with  the 
profiles  and  the  spectral  plots  as  a  group  in  order  to  understand  the  boundar>'  layer  be- 
havior. 

Notice  how  httle  the  turbulence  pulse  affects  the  spectral  distribution  ver\'  near 
the  surface  in  Figure  38.  The  velocity  spectra  for  the  undisturbed  flow  and  recovery 
period  regimes  are  ver>"  similar  and  the  total  power  distribution  curves  for  the  two  re- 
gimes are  comparable  in  shape  and  magnitude.  In  addition,  there  is  no  unusual  recovery 
period  in  this  region.  The  turbulence  pulse  comes  and  goes  without  affecting  the  spectral 
characteristics  of  the  fiow  ver\"  near  the  surface.  Higher  in  the  boundan."  layer, 
Figure  40  and  Figure  41,  the  high  frequency  components  become  more  significant  in 
the  undisturbed  fiow  and  recovery  period  regimes  while  the  total  power  contained  in 
both  regimes  drops  sharply.  This  indicates  more  spectral  power  (turbulence)  is  con- 
tained in  the  boundary  layer  very  near  the  surface  than  in  the  intermediate  regions  or 
even  in  the  external  fiow  region.  Note:  The  author  is  unable  to  explain  the  unusual 
shapes  of  the  undisturbed  fiow  velocity  spectra  in  Figure  41  and  Figure  42  or  the  ap- 
parent increase  in  undisturbed  fiow  total  power  in  Figure  42. 

Together,  the  velocity  spectra  and  total  power  distribution  plots  show  that  the 
fiow  recovers  ver\'  quickly  from  the  turbulence  pulse  and  the  pulse  does  not  significantly 
affect  the  characteristics  of  the  recovery  period.  In  other  words,  the  undisturbed  fiow 
and  recoverv'  period  velocity  spectra  and  total  power  distribution  curves  are  virtually 
identical  at  all  positions  in  the  boundary  layer. 

The  author  believes  an  understanding  of  the  mechanics  operating  in  the 
boundary  layer  at  this  chord  location  is  essential  in  understanding  the  more  complex 
responses  observed  at  50%  and  70%  chord.  The  discussion  that  follows  is  the  author's 
interpretation  of  the  operating  mechanics. 
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38.      Velocity  Spectra  and  Total  Power  Distribution,  30%  Chord,  y  =  0.0032' 
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Figure  39.      Velocity  Spectra  and  Total  Po\^er  Distribution,  30%  Chord,  y=  0.0071" 
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Figure  40.      Velocity  Spectra  and  Total  Power  Distribution,  30%  Chord,  y  =  0.0184" 
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Figure  41.      Velocity  Spectra  and  Total  Power  Distribution,  30%  Chord,  y=  0.0372' 
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Figure  42.      Velocity  Spectra  and  Total  Power  Distribution,  30%  Cliord,  y  =  0.09-14' 
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The  data  show  a  normally  laminar  flow  being  uniformly  disturbed  by  a  strong, 
highly  turbulent  pulse  from  the  external  flow.  The  boundar\'  layer  responds  to  the 
passing  pulse  by  momentarily  transitioning  to  extremely  turbulent  flow.  Once  the  pulse 
has  passed,  the  boundar\'  layer  rapidly  returns  to  low  turbulence  laminar  flow.  Note 
however,  that  the  velocity  time  history  graph  (Figure  36)  shows  that  the  mean  velocity 
quite  near  the  surface  increases  in  response  to  the  turbulence  pulse.  Furthermore,  this 
velocity  increase  remains  in  evidence  for  a  considerable  period  after  the  velocity  profile 
returns  to  laminar.   The  cause  of  this  velocity  increase  is  important. 

The  turbulence  time  histories  show  there  is  considerable  turbulence  and  mixing 
throughout  the  entire  boundary*  layer  in  response  to  the  pulse.  This  seems  to  indicate 
that  momentum  transfer  down  through  the  boundary  layer  is  possibly  responsible  for  the 
velocity  increase  near  the  surface.  The  velocity  time  histories  seem  to  support  this 
supposition  with  regard  to  observed  characteristics  of  the  velocity  deficit.  The  growth 
of  the  velocity  deficit  as  the  surface  is  approached  is  consistent  with  momentum  transfer 
down  through  the  laver. 

The  boundary  layer  region  ver\"  near  the  wing  surface  cannot  transfer  momen- 
tum due  to  the  solid  wing  boundarv^  Thus,  the  near-surface  resion  absorbs  all  the 
transferred  momentum  as  a  velocity  increase.  This  is  evidenced  by  the  sharp  decrease 
of  the  velocity  deficit  and  rapid  velocity  increase,  or  flow  acceleration,  just  following  the 
pulse.  These  two  phenomena  are  clearly  observed  in  the  bottom  two  velocity  time  his- 
tories. Finally,  the  velocity  increase  remains  in  the  near-surface  region  due  to  the  posi- 
tive velocity  gradient  ( — — )  in  the  boundar\'  layer  and  is  slowly  dissipated  by  the 
viscous  eflects  of  the  laminar  laver. 

In  summar}",  the  velocity  increase  in  the  near-surface  region  of  the  boundar\" 
layer  appears  to  be  the  result  of  a  downward  transfer  of  momentum  within  the  boundary 
layer  induced  by  the  high  turbulence  of  the  turbulence  pulse. 
2.     50%  Chord 

Figure  43  shows  the  ensemble  average  velocity  and  turbulence  profiles  at  50% 
chord.  The  boundary"  layer  response  at  this  position  is  quite  unlike  the  simple  laminar- 
turbulent-laminar  response  observed  at  30%  chord.  To  begin  with,  the  undisturbed  flow 
at  sample  150  is  transitional/turbulent  with  5%  to  6%  turbulence  near  the  surface.  The 
flow  becomes  extremely  turbulent  at  sample  350  with  a  high  level  of  turbulence  across 
the  entire  profile.  During  the  recover}-  period,  the  flow  becomes  what  appears  to  be  fully 
laminar  by  sample  650;  the  turbulence  is  extremely  low  over  the  entire  region  and  the 
boundary  layer  thickness  has  decreased  from  around  0.08  inches  at  sample  150  to  0.04 
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inches.  Samples  750  and  850  show  the  flow  returning  to  the  original 
transitional/turbulent  state  as  evidenced  by  increasing  levels  of  turbulence  in  the  near 
surface  resion. 

The  smoothed  velocity  time  histories  (Figure  44)  clearly  show  the  influence  of 
the  velocity  deficit  over  the  entire  boundary  layer.  The  corresponding  smoothed  turbu- 
lence time  histories  (Figure  45),  however,  do  not  show  any  turbulence  penetration  into 
the  near-surface  region.  Apparently  the  high  level  of  natural  near-surface  turbulence 
overshadows  the  effects  of  the  turbulence  pulse.  The  apparent  answer  to  this  dilemma 
lies  in  the  velocity  spectra  and  total  power  distribution  plots  of  the  undisturbed  flow  re- 
gime and  the  turbulence  pulse  regime  near  the  surface. 

Shown  in  Figure  46  are  the  velocity  spectra  and  total  power  distribution  curves 
taken  near  the  bottom  of  the  boundary  layer.  This  figure  corresponds  to  the  near- 
surface  region  of  the  mean  velocity  and  mean  turbulence  traces.  The  undisturbed  flow 
total  power  trace  is  quite  broad  banded  with  significant  power  contributions  from  the 
frequency  components  up  to  4  KHz  due  to  the  high  levels  of  turbulence  generated  near 
the  surface.  The  turbulence  pulse,  on  the  other  hand,  is  a  low  frequency  dominated 
phenomenon  having  no  significant  spectral  power  components  above  1  KHz.  Thus, 
while  the  turbulence  intensity  of  the  undisturbed  flow  and  turbulence  pulse  are  almost 
equal  in  magnitude,  the  spectral  makeup  of  each  regime  varies  greatly.  The  distinct  na- 
ture of  the  turbulence  pulse  makes  it  felt  throughout  the  turbulent  near-surface  layers. 
The  effect  of  the  turbulence  pulse  on  the  undisturbed  flow  is  clearly  seen  in  this  figure. 
The  undisturbed  flow  spectrum  is  altered  by  the  turbulence  pulse  such  that  the  recovery 
period  spectrum  resembles  the  turbulent  pulse  spectrum. 

Figure  47  through  Figure  51  depict  the  changes  in  the  velocity  spectra  and  total 
power  distribution  at  increasing  vertical  distances  in  the  boundan,'  layer.  Near  the  sur- 
face (Figure  46),  the  turbulence  pulse  clearly  alters  the  recovery  period  spectrum  by  re- 
moving the  high  frequency  components  contained  in  the  undisturbed  flow  regime.  In 
the  middle  of  the  boundary  layer  (Figure  49),  the  recovery  period  spectrum  begins  to 
resemble  that  of  the  undisturbed  flow  showing  the  decreasing  impact  of  the  turbulence 
pulse  on  the  flow.  Finally,  in  the  external-flow  region  (Figure  51),  the  undisturbed  flow 
and  recovery  spectra  match,  indicating  that  the  efTects  of  the  turbulence  pulse  are  purely 
transient.  Of  interest  are  the  very  large  spikes  in  the  spectra  of  the  recover}'  period  plots 
for  y<  0.1 250  inches.  The  author  speculates  that  these  are  Tollmien-Schlichting  waves 
generated  during  the  laminar  to  turbulent  transition. 
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Figure  46.      Velocity  Spectra  and  Total  Po^ver  Distribution,  50%  Chord,  y  =  0.0032' 
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Figure  47.      Velocity  Spectra  and  Total  Po»er  Distribution,  50%  Chord,  y=  0.0083' 
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Figure  48.      Velocity  Spectra  and  Total  Poner  Distribution,  50%  Chord,  y=  0.0236' 
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Figure  49.      Velocity  Spectra  and  Total  Po^^er  Distribution,  50%  Chord,  y=  0.0486' 
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Figure  50.      Velocity  Spectra  and  Total  Power  Distribution,  50%  Chord,  y=  0.1250' 
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Figure  51.      Velocity  Spectra  and  Total  Power  Distribution,  50%  Chord,  y=  0.2297' 
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The  total  power  distribution  plots  show  trends  similar  to  the  velocity  spectra 
curves.  Near  the  surface  (Figure  46),  the  turbulent  power  is  quite  high  for  all  three  flow 
regimes.  Although  the  turbulent  pulse  does  not  affect  the  recovery  period  power  level, 
the  frequency  distribution  of  the  power  is  clearly  altered  by  the  turbulence  pulse.  The 
undisturbed  flow  curve  rises  slowly  to  the  asymptote  while  the  recovery  period  curve 
rises  sharply  and  levels  off  very  quickly.  Near  the  middle  of  the  boundary  layer 
(Figure  48  and  Figure  49),  the  recovery  period  turbulence  power  is  much  lower  than 
that  of  the  free  stream.  This  shows  that  there  is  a  definite  recovery  process  operating  in 
this  region  as  the  flow  is  energized  back  to  the  free  stream  condition.  In  the  external- 
flow  region  (Figure  51),  the  undisturbed  flow  and  recovery  period  cumulative  power  are 
nearly  equal  indicating  that  the  effects  of  the  turbulence  pulse  are  purely  transitory. 

The  reverse-transition  to  laminar  flow  (variously  called  reversion,  inverse  or  re- 
verse transition,  relaminarization,  and  laminarization)  observed  in  Figure  43  is  identical 
to  the  phenomenon  observed  by  Howard  [Ref  6]  for  a  wing  immersed  in  propeller  wash. 
It  appears  that  the  turbulence  pulse  triggers  some  mechanism  in  the  boundary  layer  to 
laminarize  the  inherently  turbulent  boundary  layer  in  this  investigation  as  the  propeller 
wash  did  in  Howard's  work. 

Narasimha  and  Sreenivasan,  in  their  comprehensive  survey  of  relaminarizing 
flows,  identified  four  major  reversion  mechanisms:  reversion  by  dissipation;  reversion  in 
stablv  stratified  flows;  reversion  in  hichlv  accelerated  flows:  and  relaminarization  bv  ro- 
tation.   [Ref  16] 

Laminarization  by  sudden  acceleration  of  the  flow  in  the  near  surface  region 
appears  to  be  the  mechanism  at  work  in  this  instance.  Analysis  of  the  turbulence  pulse 
regime  of  Figure  44  reveals  instantaneous  accelerations  of  8,000  to  10,000  ft/sec^  are  felt 
by  the  flow  near  the  surface.  The  flow  appears  to  be  accelerated  both  by  the  turbulence 
pulse  itself  during  recovery  from  the  velocity  deficit  and  by  momentum  transfer  down 
through  the  boundarv"  layer  (see  previous  section). 

One  of  the  most  widely  used  parameters  that  have  proven  to  predict  the  onset 
and  course  of  laminarization  of  a  turbulent  boundar>'  layer  is: 

K  =    -^ -^  (10) 

where  K  is  the  acceleration  parameter,  Ue(x)  is  the  edge  free  stream  velocity,  v  the 

kinematic  viscositv,  and  — - — -  the  instantaneous  velocitv  derivative.     K,  as  used  in 

dx 
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Equation  10,  is  purely  a  free  stream  parameter  that  carries  no  information  at  all  about 
the  boundary'  layer  itself  Numerous  investigators  have  proposed  that  values  of  K 
greater  than  2  x  10"*  to  3.5  x  lO"'  will  give  rise  to  laminarization  in  low  speed 
incompressible  flows.  [Ref  16] 

In  this  investisation,  the  accelerated  flow  is  observed  in  the  boundarv  laver  and 
in  the  external  free  stream  as  well.  For  this  reason,  the  author  chose  to  evaluate  K  lo- 
cally within  the  boundary  layer  as  a  means  of  predicting  laminarization.  Applied  in  this 
manner,  K  becomes: 

V        dU(x,v)  ,    ^ 

L-'(x,y)        dx 

where  L'(x,y)  is  the  local  ensemble  average  velocity  at  height  y  in  the  boundary  layer. 

Values  of  the  acceleration  parameter,  K,  range  from  2.2  x  10"*  at  the  near- 
surface  region  to  1.7  x  10"*  in  the  upper  regions  of  the  boundary  layer  at  50%  chord 
using  local  velocity  values.  The  value  calculated  near  the  surface  is  clearly  sufficient  to 
induce  laminarization  of  the  turbulent  boundary  layer  at  this  position.  It  is  interesting 
to  note  that  K  is  much  higher  at  30°  o  chord  where  it  remains  fairly  constant  at  4.4  x 
10"*  over  the  entire  boundary  layer  region;  but  K  does  not  play  a  significant  role  at  30% 
chord  other  than  hastening  the  return  to  naturally  laminar  flow. 
3.     70%  Chord 

The  ensemble  average  velocity  and  turbulence  profiles  plotted  in  Figure  52 
show  yet  entirely  different  characteristics  than  observed  at  30°/o  and  50%  chord.  The 
undisturbed  flow  is  fully  turbulent  with  a  high  level  of  turbulence  spanning  the  entire 
boundary  layer  region  (sample  200).  The  turbulence  pulse  flattens  the  velocity  profile 
while  inducing  high  turbulence  (6%  -  8%)  uniformly  over  the  vertical  region  (samples 
350  and  410).  During  recover)',  the  velocity  profiles  become  laminar  like  (sample  720) 
as  the  turbulence  decreases  in  the  middle  and  upper  regions  of  the  boundary  layer.  The 
turbulence  traces  appear  transitional/turbulent  in  the  recovery  regime.  Sample  850 
shows  the  flow  becoming  fully  turbulent  again  with  the  boundar>'  layer  height  increasing 
and  with  increasing  turbulence  in  the  lower  and  middle  section  of  the  boundarv  laver. 
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Figure  52. 


Ensemble  Average  Profiles,  70%  Chord 
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Note  that  the  turbulence  level  near  the  surface  remains  high  and  fairly  constant 
throughout  the  period.  Clearly  the  laminarization  process  is  at  work  here  also,  but  the 
effects  are  greatly  reduced. 

The  smoothed  velocity  time  histories  plotted  in  Figure  53  show  that  the  velocity 
deficit  is  felt  throughout  the  boundary'  laver.  The  velocitv  deficit  grows  near  the  mid- 
layer  and  decreases  again  near  the  surface.  Of  interest  in  the  velocity  time  histories  is 
the  nature  of  the  velocity  increase  following  the  turbulence  pulse.  The  velocity  increase 
is  much  higher  at  70%  chord  than  at  the  previous  wing  locations.  Note  also  the  long 
recover}"  period  and  slow  rise  and  decay  time  associated  with  the  velocity  increase. 

Figure  54  shows  the  smoothed  turbulence  time  histories.  Notice  the  high  tur- 
bulence levels  that  persist  in  the  near-surface  regions  during  the  entire  period.  The  trace 
at  y  =  0.0702  inches  would  tend  to  indicate  that  the  turbulence  pulse  affects  the  mid- 
layer  region  more  than  the  near-surface  region  at  this  chord  location. 

The  velocity  spectra  and  total  power  distribution  graphs  are  shown  in 
Figure  55  through  Figure  60.  In  the  near-surface  region  (Figure  55  through 
Figure  57),  the  turbulence  pulse  mainly  affects  the  velocity  spectrum  of  the  recoven.' 
period  while  total  power  remains  high  and  nearly  equal  over  the  three  flow  regimes. 
Toward  the  middle  of  the  layer  (Figure  5S  and  Figure  59),  the  total  power  levels  of  the 
undisturbed  flow  and  recovery  period  regimes  drop  and  the  velocity  spectra  align  more 
closely.  It  is  interesting  to  note  that  the  total  power  of  the  recover}'  period  is  less  than 
that  for  the  undisturbed  flow  in  the  middle  boundary  layer  region.  This  indicates  that 
there  is  a  definte  recover}'  process  at  work  here.  In  the  external-flow  region 
(Figure  60),  the  turbulence  pulse  does  not  affect  the  velocity  spectrum  of  the  recovery 
period.  Here  the  power  curves  for  the  undisturbed  flow  and  recover}'  period  regimes  are 
ver}'  low  magnitude  and  equal. 

The  local  acceleration  parameter  (K)  from  Equation  11,  varies  from  1.3  x  10^* 
near  the  surface  to  1.9  x  10"^  at  the  middle  of  the  boundar}^  layer.  The  low  value  of  K 
near  the  surface  indicates  that  this  region  should  not  laminarize.  This  is  borne  out  by 
the  data;  the  near-surface  region  is  dominated  by  high  turbulence  and  the  mean  velocity 
responds  very  slowly  to  the  turbulence  pulse.  Conversely,  the  middle  boundar}'  layer 
region,  with  a  value  of  K  close  to  that  necessary  for  laminarization,  has  a  laminar-like 
response:  a  sharp  drop  in  turbulence  following  the  turbulence  pulse  with  a  corresponding 
sharp  velocity  increase. 

Apparently  the  turbulence  pulse  prompts  the  middle  boundar}'  layer  toward 
laminar  flow  but  the  pulse  is  unable  to  fully  penetrate  into  the  thick,  highly  turbulent 
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Figure  55. 


Velocity  Spectra  and  Total  Power  Distribution,  70%  Chord,  y=  0.0048' 
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Figure  56.      Velocity  Spectra  and  Total  Po>>er  Distribution,  70%  Chord,  y=  0.015-1' 


83 


9 
o 
C 

•o 
« 

•e 

3 

n 
■•6 

c 


■D 
O 

•c 
«l 

Q. 

W 

> 

o 
u 

» 


O 
O 

ro 
o 


o 
o 
o 
eo 

o 

o- 

o 

o 
,o 
o 


o 
.o 
o 

CM 


Itiiiii  I  I — [iiiiirr  I — iiiniTrT — iiiiiii  i  i    iiiiiii  i  i 


I 
o 


2««(Das/i»)    Jamod    |o^Oi 
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58.      Velocity  Spectra  and  Total  Power  Distribution,  70%  Chord,  y  =  0.0702' 
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Figure  59.      Velocity  Spectra  and  Total  Power  Distribution,  70%  Chord,  y  =  0.1941" 
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Figure  60.      Velocity  Spectra  and  Total  Pouer  Distribution,  70%  Chord,  y=  0.3-193' 
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near-surface  region.  Regardless,  the  turbulence  pulse  appears  able  to  prompt  a 
laminar-like  response  at  this  chord  location,  contrar\-  to  what  one  would  intuitively  ex- 
pect. 
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V.     CONCLUSIONS 

The  time-var>'ing  eflects  of  periodic  free  stream  turbulence  have  been  investigated 
on  a  wing  boundary"  layer  at  three  representative  chord  locations  (laminar, 
transitional  turbulent,  and  turbulent)  on  the  wing.  An  experimental  method  was  devel- 
oped and  tested  to  generate  non-thrusting,  uniform  periodic  turbulence  pulses  in  the 
external  flow  over  the  test  wing.  Although  the  turbulence  pulses  were  restricted  to  one 
test  condition,  50  pulses  per  second  and  a  turbulence  intensity  of  near  10%,  they  dem- 
onstrated the  dynamic  character  of  the  boundar>'  layer  subjected  to  single-event  turbu- 
lence.  The  results  of  this  investigation  are  detailed  below. 

The  turbulence  pulse  was  a  short  (5  millisecond),  high  intensity  (10.5%  Tu),  low 
frequency  dominated  (no  significant  spectral  components  above  1  KHz),  non-thrusting 
event. 

The  turbulence  pulse  was  shown  to  fully  penetrate  the  boundary  layer  at  all  three 
chord  locations.  It  was  observed  that  the  turbulence  pulse  suppressed  the  high  fre- 
quency spectral  components  in  the  recover}'  period  regime.  Actual  penetration  was  less 
near  the  surface  in  the  turbulent  flow  region  of  the  wing,  but  the  effects  of  the  pulse  were 
still  quite  evident. 

The  boundar\-  layer  time-response  showed  a  distinct  transition  from  undisturbed 
(laminar,  transitional,  or  turbulent)  flow  to  highly  turbulent  flow  when  subjected  to  the 
turbulence  pulse  at  all  three  chord  locations. 

The  recovery  process  leading  back  to  undisturbed  flow  was  highly  dependent  on  the 
nature  of  the  undisturbed  flow  as  shown  below: 

Laminar:  Quick  boundary  layer  response  with  a  rapid  return  to  undisturbed  laminar 
flow. 

Transitional  Turbulent:  Slower  boundan.^  layer  response  than  in  the  laminar  region. 
Laminarization  was  evident  during  the  recovery  process  followed  by  a  slow  return  to 
undisturbed  transitional  turbulent  flow.  The  boundarv'  layer  thickness  was  decreased 
during  the  recovery  period,  and  the  turbulence  intensity  dropped  to  1,6  of  its  normal 
undisturbed  value  near  the  surface. 

Turbulent:  Slow  boundarv'  layer  response.  Laminarization  is  quite  evident  in  the  mid 
and  upper  portions  of  the  boundary  layer  and  evident  to  a  lesser  degree  near  the  sur- 
face. There  was  a  marked  decrease  in  boundarv  laver  thickness  during  the  recovcrv 
period  in  this  region  before  returning  to  the  undisturbed  turbulent  value. 
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Laminarization  appears  to  be  mainly  prompted  by  rapid  acceleration  of  the  flow 
internal  to  the  boundary'  layer.  The  internal  flow  is  accelerated  both  by  the  turbulence 
pulse  itself  during  recovery  fl^om  the  velocity  deficit  and  by  momentum  transfer  down 
through  the  boundary  layer.  The  velocity  time-history  traces  show  that  the  small  ve- 
locity deficit  associated  with  the  turbulence  pulse  in  external  flow  is  amplified  in  the 
middle  boundar\-  layer  region  due  to  momentum  transfer  down  through  the  layer.  The 
near-surface  flow  experiences  a  rapid  velocity  increase  as  a  result  of  the  momentum 
transfer. 

The  total  spectral  power  content  of  the  flow  regimes  (undisturbed,  flow,  turbulence 
pulse,  and  recover^'  period)  was  more  dependent  on  vertical  position  in  the  boundary' 
layer  than  on  chord  location.  Near  the  surface,  the  total  power  of  all  three  regimes  was 
ver\'  high  and  nearly  equal  even  though  the  spectral  distribution  of  the  power  varied 
considerably.  The  middle  boundary  layer  region,  on  the  other  hand,  is  characterized  by 
a  high  power  turbulence  pulse,  medium  power  level  undisturbed  flow,  and  a  low  power 
recovery  period.  The  very  low  power  level  observed  in  the  recovery  period  shows  that 
the  turbulence  pulse  drives  the  flow  to  a  ver\'  low  turbulent  energy  state.  Laminarization 
occurs  as  there  is  insufficient  spectral  power  to  support  turbulent  flow  at  this  point.  The 
flow  slowly  gains  power  and  transitions  back  to  the  undisturbed  flow  regime  as  the 
steady  state  characteristics  once  again  prevail.  At  the  edge  of  the  boundary  layer,  the 
power  level  of  the  turbulence  pulse  remains  quite  high  while  the  power  levels  of  the  un- 
disturbed flow  and  recovery  regimes  are  low  and  nearly  equal. 
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VI.     RECOMMENDATIONS 

This  investigation  demonstrated  the  stabilizing  nature  of  periodic  turbulence  on  the 
transitional  and  turbulent  regions  of  the  wing.  The  study  was  done  at  a  single  low  angle 
of  attack  at  three  representative  chord  locations  (laminar,  transitional,  and  turbulent). 
Further  experimental  work  should  be  done  at  additional  chord  locations  and  at  higher 
angles  of  attack  to  further  study  the  stabilizing  effect  of  periodic  turbulence  on  attached 
boundary'  layers  and  on  flow  near  separation. 

More  research  also  needs  to  be  undertaken  to  better  understand  the  mechanism  at 
work  in  the  laminarization  process.  Future  studies  should  investigate  the  individual 
pulse  passage  cycles  in  greater  detail  to  determine  if  ensemble-averaging  techniques  ad- 
equately represent  the  highly  dynamic  laminarization  process,  or  if  statistical  methods 
or  a  single  pulse  passage  cycle  analysis  would  be  more  instructive. 

Data  collection  was  extremely  time  consuming  with  the  present  configuration.  The 
following  modifications  should  be  implemented  for  future  investigations: 

Store  and  process  all  data  in  binar\-  integer  format  rather  than  in  ASCII  real  format. 
This  will  reduce  the  mass  data  storage  requirement  by  a  factor  of  nine  and  also  greatly 
reduce  data  transfer  and  processing  time. 

Expand  the  random  access  memor\"  of  the  data  acquisition  computer  so  that  the  data 
collection  and  traverser  control  software  routines  can  be  managed  by  one  personal 
computer.   This  will  allow  for  full  automation  of  the  test  sequence. 

Use  a  hot-film  probe  instead  of  a  hot-wire  probe  in  the  NFS  low  speed  wind  tunnel 
to  avoid  contamination  induced  drift  problems. 

Decrease  the  digitizing  rate.  This  investigation  has  shown  that  little  useful  spectral 
information  is  contained  in  the  frequencies  above  10  KHz.  As  such,  the  sampling  rate 
could  easily  be  reduced  to  20  KHz  with  no  loss  of  meaningful  data.  This  step  alone 
will  reduce  mass  storage  by  60'^/b  with  a  comparable  reduction  in  data  handling  and 
processing  time. 

Future  research  should  also  be  directed  toward  investigating  the  effect  of  periodic 

or  single-event  turbulence  on  an  airfoil  near  stall  or  in  post-stall  flight  with  regard  to  loss 

of  attached  flow  and  loss  of  control.   Questions  to  be  answered  include: 

What  is  the  effect  of  the  observed  stabilizing  phenomena  on  flow  near  separation? 

Can  the  results  be  worked  to  an  advantage  at  high  angles  of  attack  with  regard  to 
aircraft  supermaneuverability  or  controllabihty?  This  is  important  as  augmented  lift 
from  dynamic-stall  maneuvering  will  have  limited  application  unless  the  pilot  can 
maintain  control  to  transition  quickly  from  one  maneuvering  state  to  another. 
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Directly  related  areas  of  interest  include  continued  studies  of  the  effects  of  wind 
gusts,  ship  airwake  turbulence,  missile  wake,  and  aircraft  wake  on  aircraft  stability  and 
controllability.  Future  studies  of  canard  wake/airfoil  interaction  at  high  lift  seem  perti- 
nent, especially  in  relation  to  the  continuing  supermaneuverability  studies  at  NPS. 
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APPENDIX  A.     AIRFOIL  COORDINATES 


UPPER 


x/c 


y/c 


LOWER 


x/c 


y/c 


1.000 

.0000 

1.000 

.0000 

.990 

.0008 

.990 

-.0002 

.980 

.0015 

.980 

-.0005 

.970 

.0022 

.970 

-.0008 

.960 

.0030 

.960 

-.0010 

.950 

.0039 

.950 

-.0011 

.940 

.0050 

.940 

-.0011 

.  930 

.0064 

.  930 

-. 0009 

.920 

.0080 

.920 

-.0007 

.910 

.0100 

.  910 

-.  0003 

.  900 

.0121 

.  900 

.0000 

.890 

.0142 

.890 

.0006 

.  850 

.0220 

.880 

.0010 

.800 

.0310 

.850 

.0023 

.  750 

.  0400 

.800 

.0042 

.  700 

.0483 

.750 

.0056 

.650 

.  0561 

.  700 

.0066 

.600 

.0635 

.650 

.  0071 

.550 

.0704 

.600 

.0072 

.500 

.0766 

.550 

.0067 

.450 

.0816 

.540 

.0065 

.400 

.0858 

.530 

.0063 

.  350 

.  0894 

.520 

.0061 

.  300 

.0918 

.510 

.  0059 

.  250 

.  0927 

.500 

.0056 

.240 

.  0926 

.490 

.0053 

.230 

.0923 

.480 

.0050 

.  220 

.  0919 

.470 

.0047 

.  210 

.0914 

.460 

.  0044 

.  200 

.0907 

.450 

.0040 

.  190 

.0900 

.440 

.  0036 

.  180 

.0888 

.430 

.0031 

.170 

.0876 

.420 

.  0027 

.  160 

.0862 

.410 

.0022 

.  155 

.  0854 

.400 

.0018 

.  150 

.0846 

.  390 

.0013 

.  145 

.0837 

.380 

.  0009 

.  140 

.0828 

.370 

.  0004 

.  135 

.0818 

.  360 

.0000 

.  130 

.0808 

.340 

-.0008 

.  125 

.0797 

.320 

-.0015 

.  120 

.0785 

.300 

-.0020 

.115 

.0773 

.280 

-.0025 

.  110 

.0760 

.260 

-.0029 

.  105 

.  0746 

.240 

-.0032 

.  100 

.0731 

.220 

-.  0034 
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095 

.0716 

090 

.0700 

085 

.  0684 

080 

.  0666 

075 

.0647 

070 

.  0628 

060 

.0585 

050 

.0537 

040 

.0481 

030 

.0415 

020 

.0335 

015 

.0286 

010 

.0233 

005 

.0169 

210 

-.0035 

200 

-.  0034 

180 

-.0033 

150 

-.0030 

100 

-.  0016 

080 

-.  0008 

060 

-.  0000 

050 

.0002 

040 

.  0003 

030 

.0000 

025 

-.0002 

020 

-.0008 

015 

-.0016 

010 

-.0008 

005 

.0020 
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APPENDIX  B.     DATA  SAMPLE  SPACING 


Sample 

30%  Chord 

50%  Chord 

70%  Chord 

Number 

(in) 

(in) 

(in) 

01 

0.  0032 

0.0032 

0.  0048 

02 

0.0033 

0.0034 

0.0060 

03 

0.  0038 

0.  0041 

0.0075 

04 

0.0046 

0.0051 

0.0095 

05 

0.0057 

0.0064 

0.0122 

06 

0.0071 

0.  0083 

0.0154 

07 

0.0087 

0.0106 

0.0192 

08 

0.0107 

0.  0132 

0. 0236 

09 

0.0130 

0.0163 

0.0286 

10 

0.0156 

0.0197 

0. 0342 

11 

0.  0184 

0.  0236 

0.0427 

12 

0.0217 

0.0278 

0.  0469 

13 

0.0251 

0.  0324 

0.0541 

14 

0.  0288 

0.0374 

0.0619 

15 

0. 0330 

0. 0428 

0.0702 

16 

0.0372 

0.0486 

0.0790 

17 

0.  0418 

0.0547 

0.  0884 

18 

0.0467 

0.0612 

0.0982 

19 

0.  0518 

0.0680 

0.  1086 

20 

0.  0572 

0.0751 

0.  1194 

21 

0. 0628 

0.0826 

0. 1307 

22 

0.0687 

0.0904 

0.  2425 

23 

0.0748 

0.  0986 

0.  1547 

24 

0.0811 

0.  1071 

0. 1674 

25 

0.0877 

0.  1158 

0. 1805 

26 

0.0944 

0.  1250 

0.  1941 

27 

0.  1014 

0.  1343 

0.  2080 

28 

0. 1087 

0. 1440 

0.2223 

29 

0.  1162 

0.  1538 

0.2371 

30 

0. 1238 

0. 1640 

0.2521 

31 

0.  1316 

0.  1743 

0.  2676 

32 

0.  1396 

0.  1850 

0.2833 

33 

0. 1477 

0.  1958 

0.2994 

34 

0.  1560 

0.2070 

0.3157 

35 

0.  1644 

0.2182 

0. 3324 

36 

0.  1731 

0. 2297 

0.3493 

37 

0.  1818 

0.2413 

0.3665 

38 

0.  1906 

0.  2531 

0. 3839 

39 

0.  1996 

0.2651 

0.4016 

40 

0.2087 

0.2722 

41 

0.2894 

42 

0.  3017 
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APPENDIX  C.     PROGRAM  TRAVERSE 

100  REM  *''-->'ti"'f*-^-*^'***'''  PROGRAM  TRAVERSE.  BAS  **^f-'f*'5'f*'s'fVcV?Vf*Vf******Vf***VfyfyrVfVfVfVf**iv 

110  REM 

120  REM  THIS  PROGRAM,  WRITTEN  BY  D.  KINDELSPIRE  AND  MODIFIED  BY  R.  RENOUD 

130  REM  (SEP  88)  CONTROLS  THE  VELMEX  8300  3-D  TRAVERSER.   THE  INITIAL 

140  REM  PORTION  OF  THE  PROGRAM  IS  DEDICATED  TO  OPERATOR-SELECTABLE  OPTIONS 

150  REM  (I.E.  MOTOR  VELOCITY,  MOTOR  ACCELERATION,  ETC.)  AVAILABLE  FOR 

160  REM  EACH  MOTOR.   DEFAULT  MOTOR  VALUES  FOR  EACH  OPTION  ARE  PROVIDED 

170  REM  IN  THE  PROGRAMMING  AND  THE  OPERATOR  CAN  SELECT  THE  DEFAULTS  FOR 

180  REM  ALL  MOTORS  WITH  ONE  KEYSTROKE.   THEW  PROGRAM  GIVES  THE  OPERATOR 

190  REM  A  CHOICE  OF  MANUAL  INPUT  (MEANING  ONE  MOTOR  MOVEMENT  FOR  EACH 

200  REM  INPUT  WITH  INPUT  IN  INCHES)  OR  COMPUTER-CONTROLLED  MOVEMENT 

210  REM  THROUGH  THE  BOUNDARY  LAYER.   THE  SINE  FUNCTION-BASED  ALGORITHIM 

220  REM  USED  TO  COMPUTE  THE  MOTOR  MOVEMENTS  THROUGH  THE  BOUNDARY  LAYER 

230  REM  GIVES  MOVEMENT  COMMANDS  IN  INCREMENTS  OF  THE  MINIMUM  MOTOR 

240  REM  MOVEMENT  STEP  SIZE  (0.000125  INCHES). 

250  REM 

1000  REM 

1010  REM     OPEN  THE  COM  PORT  AND  INITIALIZE  THE  MOTOR  SETTINGS 

1020  OPEN  "coml: 1200, n, 8,1, rs,cs,ds,cd"  FOR  RANDOM  AS  in 

1030  REM  SET  MOTOR  DEFAULT  VALUES 

1040  DATA  1000,1000,1000,2,2,2,0.000125,0.000125,0.000125 

1050  READ  VI,  V2,  V3 ,  Rl,  R2,  R3,  CI,  C2,  C3 

1060  DIM  Z(IOO),  Zl(lOO),  HWV(IOO),  RMS(IOO),  VEL(IOO),  NDVEL(IOO),  VRMS(IOO) 

1070  PI  =  3. 141592654# 

1080  PRINT 

1090  PRINT 

1100  PRINT 

1110  PRINT 

1120  REM 

1130  PRINT 

1140  PRINT 

1150  PRINT 

1160  PRINT 

1170  PRINT 

1180  PRINT 

1190  PRINT 

1200  PRINT 

1210  PRINT 

1220  PRINT 

1230  PRINT 

1240  PRINT 

1250  PRINT 

1260  PRINT 

1270  PRINT 

1280  PRINT 

1290  PRINT 

1300  PRINT 

1310  PRINT 

1320  PRINT 

1330  INPUT  "MANUAL  CONTROL  OR  COMPUTER  CONTROL  (ENTER  'MAN'  or  'CP')";  CON$ 


.t-»f-«.»- «,»«.»-. 


.  »•-  »v  J,  »»^  j«»t.. . 


'.  ^.  mJ^  J>»  ^.  .f.  ><<  ^.  J,  .1.^.  .««.!•  J^  **^ 


''"•-  USER  MUST  SELECT  'CAPS  LOCK'  FUNCTION  ''"'"" 


.  »t*  J-  »'-  J-  -.',  »t,  J^  «'-  *f*  fc',  -'«  .»-  .•-  *t,  »t- . 


I  %J^  %J^  ^m  tj^  mj^  fr*4  fc'4  ^*4  *J^  *Jm  < 


L«t- J- «»-«»-  t  ' 


DISPLAY  MOTOR  DEFAULT  SETTINGS 

INITIALIZED  VALUES  FOR  ALL  MOTOR  SETTINGS:" 
VELOCITY  =  1000  STEPS/SEC" 

RAMP( MOTOR  ACCEEJRATION)  =  2   (6000  STEPS/SEC 
DEFAULT  INCREMENTAL  UNITS  ARE  INCHES" 

DO  YOU  WANT  TO  INPUT  YOUR  OWN  INCREMENT  EACH  TIME  OR  HAVE  THE" 
COMPUTER  COMPUTE  MOVEMENTS  THROUGH  A  BOUNDARY  LAYER  FOR  YOU  ?" 

FOR  COMPUTED  BOUNDARY  LAYER  MOVEMENT,  YOU  MUST  KNOW:" 

1.  WHICH  MOTOR  WILL  MOVE  THE  APPROPRIATE  DIRECTION. " 

2.  HOW  THICK  (IN  INCHES)  IS  THE  B.  L.  " 

3.  HOW  MANY  DATA  COLLECTION  POINTS  DO  YOU  WANT  IN  THE  B.  L.  " 

NOTE!  !  USE  MANUAL  CONTROL  TO  INITIALIZE  PROBE  POSITION  BEFORE" 
SELECTING  COMPUTER  CONTROLLED  MOVEMENT.      " 
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1340  IF  CON$  =  "CP"  THEN  3130 

1350  REM  OPTION  TO  CHANGE  DEFAULT  SETTINGS  OF  VELOCITY  OR  ACCELERATION  RAMP 

1360  PRINT 

1370  PRINT 

1380  PRINT  "  DO  YOU  WANT  TO  CHANGE  THE  VELOCITY  OR  ACCELERATION  RAMP" 

1390  PRINT  "        DEFAULT  SETTINGS?  (Y  or  N)" 

1400  PRINT 

1410  PRINT  "IF  'NO',  THIS  PROGRAM  WILL  THEN  LET  YOU  DEFINE  THE" 

1420  PRINT  "DISTANCE  YOU  WANT  TO  MOVE  (IN  INCHES).  IF  'YES'," 

1430  PRINT  "YOU  CAN  CHANGE  ANY  OR  ALL  OF  THE  DEFAULT  SETTINGS  FOR  ANY  MOTOR." 

1440  PRINT 

1450  PRINT 

1460  PRINT 

1470  INPUT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THE  DEFAULT  SETTINGS?  (Y  or  N)";  D$ 

1480  IF  D$  =  "Y"  THEN  1520 

1490  IF  D$  =  "N"  THEN  2150 

1500  REM 

1510  REM   ***•>''   OPERATOR  SELECTED  MOTOR  VARIABLES  **'V,v* 

1520  PRINT 

1530  PRINT 

1540  INPUT  "WHICH  DEFAULT  VALUE?  (ENTER  'I'FOR  VELOC  OR  '2'  FOR  ACCEL  RAMP)"; 

1550  ON  L  GOTO  1620,  1860 

1560  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  DEFAULT  VELOCITY?  (Y  OR  N)" 

15  70  INPUT  V$ 

1580  IF  V$  =  "Y"  THEN  1620 

1590  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  DEFAULT  ACCELERATION  RAMP?  (Y  or  N)" 

1600  IF  R$  =  "Y"  THEN  1920 

1610  IF  R$  =  "N"  THEN  1380 

1620  PRINT 

1630  PRINT 

1640  INPUT  "WHICH  MOTOR  VELOCITY  DO  YOU  WISH  TO  CHANGE?  (1,2,  or  3)";  J 

1650  ON  J  GOTO  1660,  1760,  1810 

1660  PRINT 

1670  PRINT 

1680  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  #l";  VI 

1690  PRINT 

1700  PRINT 

1710  PRINT 

1720  PRINT  "DO  YOU  WANT  TO  CHANGE  VELOCITY  OF  ANOTHER  MOTOR?  (Y  OR  N)" 

1730  INPUT  V$ 

1740  IF  V$  =  "Y"  THEN  1620 

1750  IF  V$  =  "N"  THEN  1360 

1760  PRINT 

1770  PRINT 

1780  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  2";  V2 

1790  PRINT 

1800  GOTO  1710 

1810  PRINT 

1820  PRINT 

1830  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  i^3'';    V3 

1840  PRINT 

1850  GOTO  1710 

1860  PRINT 

1870  PRINT 

1880  INPUT  "WHICH  MOTOR  ACCEL  RAMP  DO  YOU  WANT  TO  CHANGE?  (1,  2,  or  3)";  K 

1890  ON  K  GOTO  1900,  1990,  2050 
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1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 


PRINT 
PRINT 

INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  /n";  Rl 
PRINT 
PRINT 

PRINT  "DO  YOU  WANT  TO  CHANGE  THE  ACCEL  RAMP  OF  ANOTHER  MOTOR?  (Y  or  N)?" 
INPUT  RM$ 

IF  RM$  =  "Y"  THEN  1860 
IF  RM$  =  "N"  THEN  1380 
PRINT 
PRINT 

INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  #2";  R2 
PRINT 
PRINT 
GOTO  1930 
PRINT 
PRINT 

INPUT  "ENTER  DESIRED.  ACCELERATION  RAMP  OF  MOTOR  y/3";  R3 
PRINT- 
PRINT 
GOTO  1930 
REM 

REM  DEFINE  DISTANCE  TO  MOVE  MOTOR 
PRINT 
PRINT 
PRINT 
REM  INITIALIZE  MOTOR  INCREMENTS  TO  ZERO 

11  =  0 

12  =  0 

13  =  0 
PRINT 

#»  #*  t\ TV  t\  *»  /»  #»  «*  *»  /»  #»  f>  #»  #t  #t  *T ««  «\ «« *v  «\  *»  **  #r *»  «v  **  tVT\  *v  #T •»7r?r *r *% 7577 *r7r *r7s  #* 7*77*57*7*  **  «*  gt  «*  «v #*  **  **7*7*  #r7* 


DEFINE  WHICH  MOTOR  YOU  WANT  TO  MOVE 

NOTE!!!   A  POSITIVE  ('+')  INCREMENT  TO  A  MOTOR 
MOVES  TRAVERSER  AWAY  FROM  THAT  PARTICULAR  MOTOR 

MOTOR  in   MOVES  THE  PROBE  UPSTREAM  AGAINST  THE  FLOW 
MOTOR  #2  MOVES  THE  PROBE  TOWARD  THE  ACCESS  WINDOW 
MOTOR  #3  MOVES  THE  PROBE  VERTICALLY  DOWNWARD 


'    ^V-VVf'5V?V>V?V*7VV-VrVr'5V?V7V*VriVVr-V7ViVVfVr?ViVVf  yVVfVriVVfVc'sVVrVfVrVr^V^VV-^V-V^V-V-VV- 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

INPUT  "WHICH  MOTOR  DO  YOU  WANT  TO  MOVE?  (1,2,  or  3)";  L 

ON  L  GOTO  2350,  2460,  2530 

PRINT 

PRINT 

PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #1?" 

PRINT  "  -V7VV-*7V7V*7VVf  (ENTER  DISTANCE  IN  INCHES)  *********" 

INPUT  II 

PRINT 

PRINT 

PRINT  "DO  YOU  WANT  TO  MOVE  ANOTHER  MOTOR  ALSO?  (Y  or  N)?" 

INPUT  C$ 

IF  C$  =  "Y"  THEN  2300 

IF  C$  =  "N"  THEN  2600 


J.J.  •  • 
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2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 


PRINT 

PRINT 

PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #2?" 

PRINT  "  --'"V,v-v-v-v-v,v,v  (ENTER  DISTANCE  IN  INCHES)  'V,vv«v,vyr,v,v*" 

INPUT  12 

PRINT 

GOTO  2410 

PRINT 

PRINT 

PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #3?" 

PRINT  "  Vov^v'«v-v,v,v*  (ENTER  DISTANCE  IN  INCHES)  ******v«v*" 

INPUT  13 

PRINT 

GOTO  2400 

PRINT 

PRINT 

REM  DISPLAY  OPERATOR  SELECTED  MOTOR  VARIABLES 

PRINT 

PRINT  "SUMMARY  OF  OPERATOR  INPUTS:  " 


PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 


MOTOR  in        VELOCITY  =  ";  VI 

ACCELERATION  RAMP  =  ";  Rl 

INCREMENTAL  DISTANCE  =  ";  II;  "INCHES" 
MOTOR  #2   VELOCITY  =  ";  V2 

ACCELERATION  RAMP  =  ";  R2 

INCREMENTAL  DISTANCE  =  ";  12;  "INCHES" 
MOTOR  #3   VELOCITY  =  ";  V3 

ACCELERATION  RAMP  =  ";  R3 

INCREMENTAL  DISTANCE  =  ";  13;  "INCHES" 


PRINT  " 
PRINT 

PRINT 

PRINT 

PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or  N)" 

PRINT 

PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.   ENTER  'Y'  TO  RETURN" 

PRINT  "TO  VARIABLE  SELECTION  SUBROUTINE.  " 

INPUT  V$ 

IF  V$  =  "Y"  THEN  1360 

GOSUB  3060 

PRINT 

PRINT 

INPUT  "DO  YOU  WANT  TO  INPUT  ANOTHER  MANUAL  MOTOR  MOVEMENT  (Y  or  N)";  M$ 

IF  M$  =  "Y"  THEN  2140 

PRINT 

PRINT  "DO  YOU  WANT  TO  INPUT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT?" 

PRINT    "  •jVVrVf'sV'sVVr^VrV-    NOTE!  !  !      *-V'V-)ViV**'!V'!V    " 

PRINT  "  ALL  PREVIOUS  MOTOR  INCREMENT  INPUTS  HAVE  BEEN  ZEROIZED.  " 

PRINT  "PROGAM  WILL  LET  YOU  CHOOSE  MANUAL  OR  CP-CONTROLLED  B. L.  MOVEMENT.  " 

PRINT  "-vvcv.-,vv.-  (IF  'NO',  THE  PROGRAM  WILL  END).  v-Vr,v-;"V" 

PRINT 

INPUT  "DO  YOU  WANT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT  (Y  or  N)";  N$ 

IF  N$  =  "Y"  THEN  1220 

PRINT 

PRINT 

PRINT 
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3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 


":C3=";  C3 


":(§" 


*■>''  NOTE  !  !  !  ** 

''"•'  COMPUTER  CONTROLLED  BOUNDARY  LAYER  ''"V 
*■>''  MOVEMENT  IS  PROGRAMMED  WITH  A  *''■ 
''"'-  DEFAULTED  NEGATIVE  MOTOR  INCREMENT  ** 
**  (i.e.  MOTOR  #3  WILL  MOVE  UPWARD  ** 
■>'■*     BY  ENTERING  A  (  +  )  B.  L.  DISTANCE).  ** 


PR  I  NT    ' '  iVVrVrVrVf  iVVriMnV  Vf  iVVf  iWfVfVf  Vr  VfiVVf  iVVf  VriVVriV* ' ' 

PRINT  "        THE  PROGRAM  HAS  ENDED. " 
PR  I  NT  "    *?V5V*?vVf!V?V'!V^v?v-)'«v-;c--v^v^v»Wf*VfV«v*-v***" 

END 

REM   *-v-'«'«vv«v  MOTOR  MOVEMENT  SUBROUTINE  *-vvov,vvr*Vc* 

PRINT  #1,  "6c":  PRINT  //I,  "E";  "C1=";  CI;  ":  C2=";  C2; 

PRINT  #1,  "Il=";  II;  ":V1=";  VI;  ":  Rl=":  Rl; 

PRINT  yn,  ":  I2=";  12;  ":  V2=";  V2;  ":R2='';  R2 

PRINT  in,  "I3=";  13;  ":V3=";  V3;  ":  R3=";  R3; 

RETURN 

PRINT 

REM  -v*,v,vvr,v,v  COMPUTER  CONTROLLED  BOUNDARY  LAYER  MOVEMENT  ******* 

PRINT 

PRINT 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT 

REM  SET  INITIAL  BL  HEIGHT  AND  NUMBER  OF  DATA  POINTS  TO  ZERO 

N  =  0 

BL  =  0 

PRINT 

INPUT  "ENTER  ESTIMATED  INITIAL  POINT  DISTANCE  ABOVE  SURFACE  (INCHES)";  ZIP 

PRINT 

INPUT  "WHAT  IS  THE  B. L.  THICKNESS  (IN  INCHES)  THAT  YOU  WANT  TO  MEASURE";  BL 

PRINT 

INPUT  "HOW  MANY  DATA  POINTS  DO  YOU  WANT  TO  SAMPLE  IN  THE  B. L.  ";  N 

PRINT 

INPUT  "WHICH  MOTOR  WILL  MOVE  YOU  THROUGH  THE  B. L.  (1,2,  or  3)";  MN 

ON  MN  GOTO  3380,  3730,  4090 

PRINT 

PRINT 

REM  MOTOR  1  CP -CONTROLLED  MOTOR  MOVEMENT 

Z(0)  =  0 

11  =  0 

12  =  0 

13  =  0 

FOR  J  =  1  TO  N 

Z(J)  =  BL  *  (1  -  SIN((PI  /  2)  *  (1  +  (J  -  1)  /  (N  -  1)))) 

STP  =  CINT((Z(J)  -  Z(J  -  1))  /  .000125)  *  .000125 

Z(J)  =  Z(J  -  1)  +  STP 

II  =  -STP 

PRINT 

PRINT 

PRINT  "    iVVr^^iViV^iV****************" 

PRINT 

PRINT  "   DATA  POINT  ";  J 

PRINT 

GOSUB  3060 

PRINT 
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3580  PRINT 

3590  IF  J  =  N  THEN  3620 

3600  INPUT  "PRESS  'ENTER'  FOR  NEXT  MOVEMENT";  MOVE$ 

3610  NEXT  J 

3620  PRINT  "ALL  MOVEMENTS  COMPLETE" 

3630  PRINT 

3640  PRINT 

3650  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER  MOVEMENT  (Y  OR  N)?" 

3660  PRINT 

3670  PRINT  "IF  'Y',  THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL  CONTROL  SUBROUTINE." 

3680  PRINT  "IF  'N',  THE  PROGRAM  WILL  END." 

3690  PRINT 

3700  INPUT  "ANOTHER  MOVEMENT";  R$ 

3710  IF  R$  =  "Y"  THEN  1290 

3720  IF  R$  =  "N"  THEN  2990 

3730  PRINT 

3740  PRINT 

3750  REM  CP -CONTROLLED  MOTOR  MOVEMENT  FOR  MOTOR  2 

3760  Z(0)  =  0 

3770  II  =  0 

3780  12  =  0 

3790  13  =  0 

3800  FOR  J  =  1  TO  N 

3810  Z(J)  =  BL  *   (1  -  SIN((PI  /  2)  *   (1  +  (J  -  1)  /  (N  -  1)))) 

3820  STP  =  CINT((Z(J)  -  Z(J  -  1))  /  .000125)  *  .000125 

3830  Z(J)  =  Z(J  -  1)  +  STP 

3840  12  =  -STP 

3850  PRINT 

3860  PRINT 

3870  PRINT 

3880  PRINT  "   Vr**v«v?v»v*y«v*******Vc******" 

3890  PRINT 

3900  PRINT  "   DATA  POINT  ";  J 

3910  PRINT 

3920  GOSUB  3060 

3930  PRINT 

3940  PRINT 

3950  IF  J  =  N  THEN  3980 

3960  INPUT  "PRESS  ENTER  FOR  NEXT  MOVEMENT";  MOVE$ 

3970  NEXT  J 

3980  PRINT  "ALL  MOVEMENTS  COMPLETED" 

3990  PRINT 

4000  PRINT 

4010  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER  MOVEMENT  (Y  OR  N)?";  R$ 

4020  PRINT 

4030  PRINT  "IF  'Y',  THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL  CONTROL  SUBROUTINE." 

4040  PRINT  "IF  'N',  THE  PROGRAM  WILL  END." 

4050  PRINT 

4060  INPUT  "ANOTHER  MOVEMENT";  R$ 

4070  IF  R$  =  "Y"  THEN  1290 


4080 

IF  R$  = 

=  "N" 

THEN 

2990 

4090 

PRINT 

4100 

PRINT 

4110 

REM  CP- 

■CONTROLLED 

MOTOR 

MOVEMENT 

FOR 

MOTOR 

3 

4120 

Z(0)  = 

0 

4130 

11  =  0 
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4140  12  =  0 

4150  13  =  0 

4160  FOR  J  =  1  TO  N 

4170  Z(J)  =  BL  *   (1  -  SIN((PI  /  2)  *  (1  +  (J  -  1)  /  (N  -  1)))) 

4180  STP  =  CINT((Z(J)  -  Z(J  -  1))  /  .000125)  '>  .000125 

4190  Z(J)  =  Z(J  -  1)  +  STP 

4200  13  =  -STP 

4210  PRI^r^ 

4220  PRINT 
4230  PRINT 
4240  PRINT  "   7V**^v*yc>v*'>v**Vc*»V'jV*******'>v*" 

4250  PRINT 

4260  PRINT  "   DATA  POINT  ";  J 

4270  PRINT 

4280  GOSUB  3060 

4290  PRINT 

4300  PRINT 

4310  IF  J  =  N  THEN  4340 

4320  INPUT  "PRESS  'ENTER'  FOR  NEXT  MOTOR  MOVEMENT";  MOVE$ 

4330  NEXT  J 

4340  PRINT  "ALL  MOVEMENTS  COMPLETED" 

4350  PRINT 

4360  PRINT 

4370  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER  MOVEMENT  (Y  OR  N)?";  R$ 

4380  PRINT 

4390  PRINT  "IF  'Y',  THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL  CONTROL  SUBROUTINE." 

4400  PRINT  "IF  'N',  THE  PROGRAM  WILL  END."  , 

4410  PRINT  ' 

4420  INPUT  "ANOTHER  MOVEMENT";  R$ 

4430  IF  R$  =  "Y"  THEN  1290 

4440  IF  R$  =  "N"  THEN  2990 

4450  END  I 
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APPENDIX  D.     PROGRAM  HWCAL 


100 

REM 

110 

REM 

120 

REM 

130 

REM 

140 

REM 

150 

REM 

160 

REM 

170 

REM 

180 

REM 

190 

REM 

200 

REM 

210 

REM 

220 

REM 

230 

REM 

VfyfVf*v«vVfyf*yr*y«vv?VfVfV«vvf  PROGRAM  HWCAL.  BAS  *-^f^c*-it*-^***********-if*****^!^t***'ifif 

This  program  is  used  to  calibrate  the  hot  wire  in  the  low  speed  wind 
tunnel.   The  hotwire  is  calibrated  against  a  pitot/static  probe 
located  in  the  free  stream  adjacent  to  the  hotwire.   Pitot/static 
output  (delta  P)  is  in  cm  of  water.   Hotwire  output  from  the  IFA-100 
is  digitally  read  by  a  Metrabyte  16f  A  to  D  board  controlled  by 
Labtech  Notebook  software.   Suggested  A  to  D  speed  is  1000  Hz  with  a 
period  of  1.  01  sec.   The  first  1000  data  points  are  used  by  this 
program.   The  program  output  file  is  directly  compatable  with 
'Grapher"  for  plotting  calibration  curves  and  determining 
conversion  equations. 


»\    «ZlfS  4\  ~\    «r7v7*    4\    4>    4\    4\    4S    4\    4^  4\    *\    *\    4\    4S    4\  ~\i%    4\i\    4\    4\    4\    4\    4\    «*    4\    4>    4*   «*.   «%    4Sty    4\    4\    t\    4% 


yiMl)  IV  Fj  rl  '^<«    '«    '«    4i^Z  4%    4\    4\    4\    4\    4%    4\    4\    4\    4*    4\  77  4\    4\    4\    4%    4S    4%    4i  4\    4\    4\    «k7«    4\    4Z  4^  4^    4^  4i  4\    4*    4^  4*    4^  4\    4\    4K    4\    4iT*  7V    ««    4\    4^  4\    4\    4\    4^  4%    4%    4*    4%    4^4*  #V<V    4%    4\    4\    ««    4* 

250  REM 

260  REM  VARIABLES: 


dp 
hwv 
vfps 
ofi 


Pitot/static  probe  delta  P  (cm  water) 
Hotwire  voltage  (volts) 
Tunnel  velocity  (feet  per  second) 
Output  file  name 


1080  INPUT 
1090  PRINT 


270  REM 

280  REM 

290  REM 

300  REM 

310  REM 

^  y  11      Kr.rl     '*  '*  **  '*  **  '*  **  '*  '*  '*  '*  '*  **  7*  **  #*  *%  «« *v «« /«««#«««  ^\  ««<*««#«#««*  7%  «r  **  *v  7*7*7*  *»  ••  **  7%  *r7c7s  7v  #*  7s  7*7*  #*  *v  *r7»7r#»  *%  *r  #*  /«  *>«#»  #*  *•  **  *•  •r 

1000  DIM  DP(20),  HWV(20) 

1010  R  =  53.3:  REM  gas  constant  for  air 

1020  G  =  32. 174 

1030  CLS  :  PRINT  "Enter  output  file  name":  INPUT  ;  OFI$ 

1040  OPEN  OFI$  FOR  OUTPUT  AS  #2 

1050  PRINT  #2,  "'cm  water'",  "'hwv'",  " '  vf^^"",' ,"%lfip5  .5'" 

1060  1=1 

1070  PRINT  :  PRINT  "Enter  cm  of  water  for  calibration  point  ";  I 

DP(I) 

HWV(I)  =  0 
1100  REM  Execute  A  to  D  conversion  of  hotwire  using  Labtech  Notebook  "go" 
1110  SHELL  "  nb  go" 

1120  REM  Labtech  Notebook  writes  output  file  to  data. prn 
1130  OPEN  "c:  nb  data  data. prn"  FOR  INPUT  AS  #1 
1140  FOR  J  =  1  TO  1000 
1150  INPUT  in,    HWV 
1160  HWV(I)  =  HWV(I)  +  HWV 
1170  NEXT  J 
1180  CLOSE  in 

1190  HWV(I)  =  HWV(I)  /  1000 

1200  CLS  :  PRINT  "hvw  for  point  ";  I;  "  =  ";  HWV(I) 
1210  INPUT  "another  calibration  point  (y/n)";  C$ 
1220  IF  C$  =  "n"  THEN  GOTO  1250 
1230  1=1+1 
1240  GOTO  1070 

1250  CLS  :  INPUT  "Enter  ambient  temperature  (Deg  F)";  T 
1260  T  =  T  +  460 
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1270  INPUT  "Enter  barometric  pressure  (in.  Hg)";  P 

1280  PRINT  "Wait  -  Writing  Output  Data  File" 

1290  P  =  P  /  29.92  *    14.7  •>''    144 

1300  RO  =  P  /  (R  ''^  T  ■>''   G) 

1310  FOR  J  =  1  TO  I 

1320  VFPS  =  SQR(2  ■■''  DP(J)  *  2.048  /  RO) 

1330  PRINT  #2,  DP(J),  HWV(J),  VFPS,  HWV(J)  2,  SQR(VFPS) 

1340  NEXT  J 

1350  CLOSE  #2 

1360  CLS  :  INPUT  "Another  calibration  run  (y/n)";  C$ 

1370  IF  C$  =  "y"  THEN  GOTO  1000 

1380  PRINT  "BYE" 

1390  END 
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APPENDIX  E.     PROGRAM  BLDATA 


100 
110 

111 

112 

113 

140 

150 

160 

170 

180 

190 

200 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1125 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 


-ki-ic-irifJficMcie-idfiti'fJc    PROGRAM    BLDATA.  BAS    *****''fif^f**********-iricic*:^riri-*ir*iri-irit 

This  driver  program  controls  the  analog-to-digital  data  acquisition 
process  using  the  software  package  LABTECH  NOTEBOOK  in  conjunction 
with  the  Metrabyte  16f  A  to  D  card. 

This  progarm  allowes  the  user  to  enter  the  number  of  data  points 
for  the  run.   It  initiates  the  data  collection  by  running  LABTECH 
NOTEBOOK  program  "go"  directly  from  the  nb  directory  via  dos 
command.   This  control  program  also  generates  incremental  names  for 
the  data  files. 


REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 

CLS 

INPUT  "Enter  number  of  data  points  for  this  wing  position  (1-99)";  N 

IF  (N  <  1)  OR  (N  >  99)  GOTO  1010 

REM  """  generating  string  segments  for  data  file  names 

B$  =  MID$(STR$(1),  2):  REM  string  number  "l" 

E$  =  MID$(STR$(N),  2):  REM  ending  string  number  "n" 

X$  =  "XXXXXX" 

EX$  =  ".DAT" 

CLS 

PRINT  "Data  files  will  be  incremented  from:" 

PRINT  (X$  +  B$  +  EX$);  "  To  ";  (X$  +  E$  +  EX$) 

PRINT 

PRINT  "Enter  common  data  file  id  (XXXXXX)  for  this  run 

PRINT  "(Six  characters  maximum  -  no  extension" 

INPUT  F$ 

IF  LEN(F$)  >  6  OR  LEN(F$)  <  1  GOTO  1120 

CLS 

PRINT  N;  "Data  files  will  be  generated  and  incremented  as  follows:" 

PRINT  (F$  +  B$  +  EX$);  "   To  ";  (F$  +  E$  +  EX$) 

PRINT 

INPUT  "Are  the  number  of  points  and  file  names  ok  (y  n)";  C$ 

IF  C$  =  "n"  GOTO  1000 

IF  C$  =  "y"  GOTO  1230 

GOTO  1190 

CLS 

FOR  I  =  1-  TO  N 

PRINT  :  PRINT  "ready  for  point  ";  I;  "  of  ";  N;  "  points":  PRINT 

PRINT  "Insert  New  Data  Disk  in  Drive  A" 

PRINT  "strike  any  key  when  ready" 

BEEP 

C$  =  INKEY$:  IF  C$  =  ""  GOTO  1290 

REM  ***  generating  incremented  data  file  name 

IF  (I  >  10)  OR  (I  =  10)  THEN  1$  =  MID$(STR$( I) ,  2) 

IF  (I  <  10)  THEN  1$  =  (MID$(STR$(0),  2)  +  MID$(STR$( I) ,  2)) 

FI$  =  (F$  +  1$  +  EX$) 

PRINT  "Doing  A  to  D  conversion  and  writing  data  file  ";  FI$ 

SHELL  "cd  nb":  REM  DOS  command  to  ensure  system  in  nb  subdirectory 

SHELL  "go":  REM  executing  labtech  notebook  via  DOS  from  NB  subdir 

Rl$  =  "ren  A:  DATA. PRN  " 
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1380  R2$  =  Rl$  +  FI$ 

1390  SHELL  R2$ 

1400  NEXT  I 

1410  INPUT  "Another  run";  C$ 

1420  IF  C$  =  "y"  GOTO  1000 

1430  PRINT  "BYE" 

1440  END 
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APPENDIX  F.     PROGRAM  CONVERT 

PROGRAM  CONVERT 

C 

C  THIS   PROGRAM  CONVERTS  THE  RAW  HOT-WIRE  VOLTAGE  AND  POSITION  SENSOR 

C  VOLTAGE   INTO  VELOCITY  DATA  AN  A  SYNCHRONIZATION  PULSE.       INPUT  DATA 

C  IS   IN  ASCII    INTEGER  FORM  WITH  THE   INTEGER  VALUES   0  TO  4095   REP- 

C  RESENTING  0  TO  5  VOLTS.      THE   BAR  SENSOR  VOLTAGE  WAS   PROCESSED  TO 

C  RETURN  A     VALUE  OF   2  AT  THE  EXACT  CENTER  OF  THE   POSTION  SENSOR 

C  PULSE   AND  A  0   ELSEWHERE.       THE   HOT-WIRE   VOLTAGE   TO  VELOCITY  CONVER- 

C  SION  CONSTANTS  ARE   OBTAINED  FROM  THE  HOT-WIRE  CALIBRATION  DATA 

C  FROM  PROGRAM  HWCAL. BAS. 

C 

I^7r#»7*  *r*»  *i*fc  ****  ^v  *»  **  •**»  #»  *»  #S  ***»  #v  0\  *s#C  0\t\0%t\  #*T*7»  #r7*  7*  *%  *r*»  ««  7*  77/S  •»  7*7^*^7*7*  /»  #»  *\  #*7*  **  *\  «»  7s  ^fTTTT?*  »%  7«  **  7»  #»  **  **  **  **  **  ** 

C 

REAL  HW^V,    PROP,    A,    B,    PPROP,    TIC 

CHARACTER^nZ   INFILL,    OUTFILE,    CH0ICE*2 

INTEGER   I,J,K 

DIMENSION  DAT  (50000,2) 

A  =  0 

B  =  0 
C 
5  PRINT    '(/A   )','    ENTER  THE   INPUT  DATA  FILE  NAME: ' 

READ   (",'(A)')    INFILE 
C 

PRINT   '(/A   )',    •    ENTER  THE  OUTPUT  DATA  FILE  NAME: ' 

READ   (-'^'(A)')    OUTFILE 
C 

OPEN   (10,    FILE=INFILE,    STATUS  =    'OLD') 

OPEN   (11,    FILE=OUTFILE,    STATUS   =    'NEW') 
C 

PRINT    '(/A   )' ,    '    CONVERSION  FROM  HWV  TO  VELOCITY   IS  OF  THE  FORM: ' 

PRINT    '(/A   )',    '    VEL  =   (A  *  HWV'"^2  +  B)   *''■   2' 

PRINT    '(/A/)' ,      '    .    .    .     ' 

WRITE   (''-,''0    '    A  =    '  ,    A 

WRITE   (",-)    '    B  =    ' ,    B 

PRINT    '(/A/)',    '    CHANGE  A  OR  B   (Y/N)    ?    ' 

READ   (",'(A)')    CHOICE 

IF   (CHOICE   .EQ.     'N')   GOTO   7 

PRINT    '(/A   )' ,    '    A  =    ' 

READ   (-,-)   A 

PRINT    '(/A   )' ,    '    B  =    ' 

READ   (*,^0    B 
C 
7  PPROP=0 

PRINT   '(/A/)',    '    READING  DATA  INTO  ARRAY   -  WAIT' 

DO   10  J  =   1,    50000 

READ   (10,   *,    END=50,    ERR=50)   HWV,    PROP 

DAT  (J,l)   =  HWV 

DAT   (J, 2)   =  PROP 
10  CONTINUE 
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c 

50    CONTINUE 

PRINT  '(/A/)',  '  DOING  DATA  CONVERSION  -  WAIT  SOME  MORE' 

DO  20  I  =  1,  J-1 

DAT  (1,1)  =  (A  *  (DAT  (1,1)  *  5.  /  4095.)  ''">  2  +  B)  *•>'■  2 

PROP  =  DAT  (1,2)  *  5.  /  4095. 

TIC  =  0. 

IF  (PROP  .GT.  1.5  .AND.  PROP  .  LT.  2.  .AND.  PROP  .  LT.  PPROP) 
+TIC  =  2. 

PPROP  =  PROP 

DAT  (1,2)  =  TIC 
20    CONTINUE 
C 

PRINT  '(/A/)',  '  WRITING  OUTPUT  FILE  -  ALMOST  DONE!' 

DO  30  I  =  1,  J-1 

WRITE  (11,80)  DAT  (1,1),  DAT  (1,2) 
80    FORMAT  (F11.4,  F5.  0) 
30    CONTINUE 
C 

WRITE  (''',*)    '  #  OF  CONVERSIONS  IS  :  '  ,  (J-1)  *   2 

CLOSE  (10) 

CLOSE  (11) 

WRITE  ('^'>)  '  ANOTHER  RUN  (Y/N)  ?' 

READ  (''S'(A)')  CHOICE 

IF  (CHOICE  .EQ.  'Y')  GOTO  5 
C 

END 
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